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Superconductivity is a phenomenon where the macroscopic quantum coherence appears 
due to the pairing of electrons. This offers a fascinating arena to study the physics of 
symmetry breaking, i.e., broken gauge symmetry. However, the important symmetries in 
superconductors are not only the gauge invariance. Especially, the symmetry properties 
of the pairing, i.e., the parity and spin-singlet / spin-triplet, determine the physical 
properties of the superconducting state. Recently it has been recognized that there is 
the important third symmetry of the pair amplitude, i.e., even or odd parity with respect 
to the frequency. The conventional uniform superconducting states correspond to the 
even-frequency pairing, but the recent finding is that the odd-frequency pair amplitude 
arises in the spatially non-uniform situation quite ubiquitously. Especially, this is the 
case in the Andreev bound state (ABS) appearing at the surface/interface of the sample. 
The other important recent development is on the nontrivial topological aspects of 
superconductors. As the band insulators arc classified by topological indices into (i) 
conventional insulator, (ii) quantum Hall insulator, and (iii) topological insulator, also 
are the gapped superconductors. The influence of the nontrivial topology of the bulk 
states appears as the edge or surface of the sample, i.e., bulk-edge correspondence. 
In the superconductors, this leads to the formation of zero energy ABS (ZEABS). 
Therefore, the ABSs at the surface/interface of the superconductors are the place where 
the symmetry and topology meet each other which offer the stage of rich physics. In 
this review, we discuss the physics of ABS from the viewpoint of the odd-frequency 
pairing, the topological bulk-edge correspondence, and the interplay of these two issues. 
It is described how the symmetry of the pairing and topological indices determine the 
absence/presence of the ZEABS, its energy dispersion, and properties as the Majorana 
fermions. Various related issues such as the Helium 3, transport of Majorana fermions, 
and tunneling spectroscopies are also briefly discussed. 

KEYWORDS: Odd-frequency pairing, Andreev bound state, Bulk-edge correspondence, Majorana fermion, 
Topological index 

1. Introduction 

Two of the most important principles in physics are symmetry and topology, both 
of which are related to the quantum numbers and energy degeneracy of states. For 
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the compact symmetry group such as the 27r rotation, the discrete quantum numbers 
result. However, even a shght breaking of the symmetry destroys these features, i.e., the 
conservation law is violated and the degeneracy is lifted. In the many-body systems, it 
often happens that the symmetry is spontaneously broken as a collective phenomenon. 
Ferromagnetism and superconductivity are two representative examples of this sponta- 
neous symmetry breaking characterized by the order parameters. The magnetization, 
which breaks the rotational symmetry in the spin space, and the pair amplitude, which 
breaks the gauge symmetry, are the corresponding order parameters, respectively. In 
sharp contrast to this, the topological quantum numbers are more robust against the 
perturbations. For example, the quantization of the magnetic flux captured by one vor- 
tex in a superconductor is quantized to be 0o — hc/{2e) even in the disordered case 
without the cylindrical symmetry. This is guaranteed by the single-valued property of 
the superconducting order parameter and the stability of its winding number. It is ev- 
ident that the winding number around a singular point ( the center of the vortex in 
the above example) does not change for the continuous deformation as long as the path 
does not cross the singular point. Therefore, the topological properties often lead to 
the robustness of the system against the weak and moderate perturbations protected 
by the topological quantum number. 

It is not always clear how to separate the two origins of the quantum numbers, 
i.e., symmetry and topology. Even more interesting is the interplay between these two 
principles. It often happens that they are entangled to give a rich variety of physical 
phenomena. In this review article, we describe how these two principles are merged into 
the novel properties of superconductors. 

1.1 Symmetry classifications of hulk superconductors and relevant quantum phenomena 
Symmetry of Cooper pair is a central topic in the physics of superconductivity. 
Here we discuss uniform bulk systems. It is established that Cooper pair is formed 
between two electrons. In accordance with the Fermi-Dirac statistics, pair amplitude, 
which is a wave function of Cooper pair, must have a sign change with the exchange of 
two electrons. Symmetry of Cooper pair is customary classified into spin-singlet even- 
parity and spin-triplet odd-parity where even (odd) refers to the orbital part of the pair 
amplitude. For example, s-wavc and d-wavc pairings belong to the former case while 
p-wave pairing belongs to the latter. With broken inversion symmetry, spatial parity is 
no more a good quantum number. In that case, spin-singlet even-parity and spin-triplet 
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odd-parity pairings can mix with each other. The mixed parity pairing has become an 
important issue in non-centrosymmetric superconductors. ^'^^ It is noted that both in 
spin-singlet even-parity and spin-triplet odd-parity pairings, the pair amplitude does 
not have a sign change with the exchange of two time variables for the two electrons 
forming the Cooper pair. However, there is a possibility that pair amplitude has a sign 
change by this operation. The latter type of pairing is so called odd-frequency pairing 
originally proposed by Berezinskii.^^ There are two possibilities of odd-frequency pairing, 
i.e., spin-singlet odd-parity and spin-triplet even-parity pairings. Thus, in the presence 
of both the inversion and spin-rotational symmetries. Cooper pairs are classified into 
(i) even-frequency spin-singlet even-parity (ESE), (ii) even- frequency spin-triplet odd- 
parity (ETO), (iii) odd- frequency spin-triplet even-parity (OTE), and (iv) odd- frequency 
spin-singlet odd-parity (OSO). In strongly correlated systems, due to the retardation 
effect of the electron interaction, there is a possibility for odd-frcqucncy pairing in the 
bulk. After the prediction by Berezinskii,^^ there have been several theoretical proposals 
about odd-frequency pairing in the bulk.^"^^'' However, up to now, odd-frequency bulk 
superconductor has not been established experimentally yet. Although it is not easy 
to realize odd-frequency gap function in uniform bulk systems, it is more promising to 
induce odd-frequency pair amplitude in the non-uniform system with lower symmetry. 
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Table I. The symmetry of pair amplitude with respect to the exchange of spins, spatial coordinates, 
and time variables for possible four classes. ESE, ETO, OTE and OSO denote even-frequency spin- 
singlet even-parity, even-frequency spin-triplet odd-parity, odd-frequency spin-triplet even-parity, and 
odd-frequency spin-singlet odd-parity. 

Bergeret, Volkov and Efetov have clarified that in ferromagnet/superconductor het- 
erostructures with inhomogeneous magnetization, odd-frequency pairing is generated 
in ferromagnet.^^-* They have predicted the long range proximity effect by the odd- 
frequency pairing. Furthermore, one of the authors (Y.T.) has revealed that the odd- 
frequency pairing is possible in inhomogeneous superconducting systems even without 
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magnetic ordering. We will show ubiquitous presence of odd-frequency pairing in 
§2.20.22) 

Symmetry of the Cooper pair influences tunneling phenomena in superconducting 
heterostructures. It is known that Andrccv bound state (ABS) is generally generated 
at the surface of anisotropic superconductor where pair potential (energy gap function) 
has a sign change on the Fermi surface. 2^"'^^) By properly taking into account of the An- 
dreev reflection in anisotropic superconductors, ^^-34) been proven that the ABS 

manifests itself as a zero bias conductance peak (ZBCP) of quasiparticle tunneling spec- 
troscopy.2^'28'35) The presence of ZBCP has been observed in tunneling spectroscopy of 
high Tc cuprates.^^"^^) It has been estabhshed that a dispersionless zero energy ABS 
(ZEABS)2^) is generated at the surface for spin-singlet d^y-waye superconductors. The 
existence of the ABS influences seriously on charge transport, ''^"^^^ spin transport 
and magnetic responses. ^^"^'^^ 

ABS is also expected for spin-triplet p-wave superconductors. It induces ZBCP in 
the tunneling spectroscopy.^^' However, proximity effect via ABS is completely dif- 
ferent between spin-singlet superconductors and spin-triplet superconductors. The ABS 
cannot penetrate into a diffusive normal metal (DN) attached to the superconductor 
in the former case, while it is possible in the latter one. The underlying physics can 
be expressed by the symmetry of the induced odd- frequency pairing. We will discuss 
these exotic phenomena in §2. 



1.2 Topological orders and topological quantum numbers 

Up to now, we have discussed the classifications of the superconductors by symmetry 
as an example of quantum states. Then, is there any other way to characterize the 
quantum phases of electronic systems ? The answer is yes, and the topological orders 
and topological quantum numbers serve this purpose. Historically, the first topological 
ordered state in condensed matter physics is the quantum Hall system, i.e., the two- 
dimensional electronic systems under strong magnetic field at low temperature, where 
the Hall conductance an is quantized to be the integer times the unit e'^/h (integer 
quantum Hall effect (IQHE)) or its fractions (fractional quantum Hall effect (FQHE)) 
(e: unit charge, h: Planck constant). Here, the topological quantum number Ci (integer, 
see eq.(48)) characterizing the topological order is related to the physical observable. 
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i.e., the Hall conductance an as 

C7H - -jC„ (1) 

which explains why an is quantized as the integer multiple of e^//i.^^^ 

Another important issue related to the IQHE is the one-dimensional channels along 
the edge of the sample. A classical picture for these edge channels is that the cy- 
clotron motion of the electrons in the bulk is reflected at the boundary of the sample, 
forming the one-dimensional motion in one-direction. This property of " one-direction" 
is called "chiral". In the quantum mechanical picture, the two-dimensional electronic 
states under magnetic field form the Landau levels, whose energies are pushed up near 
the edges of the sample by the confining potential. As a consequence, some of the edge 
channels cross the Fermi energy, which carry current. A remarkable fact is that the 
direction of the propagating waves along these edge channels is one-way, i.e., chiral, 
and those with the opposite direction are on the other side of the sample, i.e., far away 
separated by the bulk. (Consider the cylindrical sample where the two edges are sepa- 
rated by the bulk.) Therefore, the backward scattering is strictly forbidden even if the 
impurity potential is there. This leads to the perfect one-dimensional conductance of 
each channel, and consequently the Hall conductance is given by 

cth = -j^n (2) 

where n is the number of the edge channels including the sign representing the chirality. 
Comparing eqs.(l) and (2), one concludes that there is a relation between the bulk 
topological quantum number (Ci) and the edge channels (n).^^'^^) This is called "bulk- 
edge correspondence" and will be explored more in details later. One can say that 
the one-dimensional electronic states are "split into halves" at the two edges of the 
sample according to the chiralities, and the fractionahzed electrons become robust. This 
fractionahzation of the electrons at the edge is realized/supported by the topological 
order in the bulk state. 

FQHE is a more complex phenomenon where the electron-electron interaction plays 
the essential role. Therefore, one needs to study the many-body wave functions, and 
the topological quantum numbers characterizing the quantum states are much more 
intriguing. The readers are referred to the textbooks for the details. ^'^'^^ However, con- 
ceptually the story is quite analogous and similar to the IQHE. The topological order 
in the bulk refiects itself as the edge channels at the boundary of the sample, i.e., the 
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bulk-edge correspondence occurs. 

A recent breakthrough is the discovery of the topological insulators which preserve 
the time- reversal symmetry T, i.e., without the external magnetic field or the sponta- 
neous magnetization.^^' With T-symmetry, the Hall response and hence the Chern 
number is zero. Therefore, it is required to define new topological quantum numbers to 
characterize the topological orders. 

The concept of topologically nontrivial states of the single-particle Hamiltonian 
can be generalized to the superconducting state as long as one considers the BCS 
Hamiltonian in the Nambu representation; 



where ijjk = (c„fct, Cnki, (^nkt^ '^nfcj,)^- ^ naive way is to consider the matrix elements 
of h{k) as that of the Bloch states discussed above. However, the situation is not so 
trivial since the degrees of freedom are redundant in ipk and ip\.. Consequently, the edge 
channels and surface states originating from the bulk-edge correspondence has the same 
redundancy, which sometimes leads to Majorana fermion. The present edge channels 
and surface states appear as the ABSs characterized by the symmetry of Cooper pair 
as discussed in the former subsection. Therefore, the physics is much richer in this 
topological superconductor than the non-interacting electron systems. 

In this review, we discuss the recent developments on the topological supercon- 
ductors as the merging point of the two principles, i.e., symmetry and topology. The 
in-depth study of this issue will clarify the interplay between these two principles in 
physics. The plan of this review follows. In §2, the odd-frequency superconducting pair- 
ing induced by the extrinsic symmetry breaking such as spatial non- uniformity and/or 
the time-reversal symmetry breaking is studied to demonstrate how the superconductiv- 
ity is infiuenced by the symmetries. In §3, we discuss the topology of the bulk states in 
superconductors and its relation to the surface/edge states (bulk-edge correspondence). 
Various related issues are described in §4. §4 is also devoted to the summary and future 
perspectives. Hereafter, we choose the unit h — ks — 1 unless explicitly written. Two 
key concepts in this review are odd- frequency pairing and edge states. In Fig.l, we 
show the relations between each subsection and above two key concepts cxphcitly. The 
readers who want to know the essence of the review urgently would be appreciated to 
look at Table IX in §4. 




(3) 



k 
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Fig. 1. The relations between eacli subsection and two key concepts, i.e., odd- frequency pairing and 
edge states, in the present review. 



2. Extrinsic symmetry breaking and odd- frequency pairing in superconduc- 
tors 

2.1 What is extrinsic symmetry breaking ? 

In this section, we discuss what is an extrinsic symmetry breaking. Typical exam- 
ples are breaking of the translational invariance or that of spin-rotational symmetry. In 
general, Cooper pair is formed between position ri and r2- Pair amplitude is a function 
oi r — (ri + r2)/2 and f — (vi — r2). When we focus on the symmetry with respect 
to the exchange of two coordinates ri and r2, i.e., f — > — f , we can define the parity 
of the Cooper pair. We employ the Fourier transformation from f to k. If we consider 
inhomogeneous superconducting systems, like vortex and junctions, translational sym- 
metry is broken. In that system, parity of the Cooper pair is no more a good quantum 
number. Then, near the vortex core or the interface / surface of the superconductor, 
mixed parity state can be realized. To be more specific, let us consider normal metal / 
superconductor (N/S) junction where the symmetry of the superconductor belongs to 
ESE. Near the interface or in the normal metal, due to the breaking of translational 
symmetry, odd-parity state can be mixed. Since the spin rotational symmetry is not 
broken, pairing symmetry of spin should be singlet. To be consistent with Fermi Dirac 
statistics, the resulting Cooper pair should be odd in frequency. Then, we can expect 
OSO paring is generated^^^ near the interface and also inside normal metal. On the 
other hand, if symmetry in the bulk is ETO, then the induced pairing symmetry by 
breaking of the translational invariance or the spatial inversion symmetry is OTE. 

Next, we discuss the breaking of the spin rotational symmetry only. By this sym- 
metry breaking, spin-singlet and spin-triplet parings are mixed with each other. Let us 
consider superconductor with ESE symmetry. In the presence of Zeeman magnetic field. 
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spin-triplet pairing state is generated. Since spatial parity is not broken, the induced 
pair amplitude should be OTE. On the other hand, if the pairing in the bulk is ETO, 
then the induced pair amplitude by breaking of spin rotational symmetry is OSO. 

From this viewpoint, paring symmetry in ferromagnet / superconductor (F/S) junc- 
tion is very interesting. Here, we specify that the symmetry of the superconductor is 
spin-singlet s-wave one which belongs to ESE. In the light of above discussion, OSO 
pairing is induced in F by translational symmetry breaking. Also by the spin rotational 
symmetry breaking, OTE (ETO) pairing is generated from ESE (OSO) pairing. Thus, 
in F, all four kinds of pairings are basically possible. In the standard case accessible 
in experiments, F is in the diffusive regime. Then, basically only s-wave pairing is pos- 
sible since it is robust against impurity scattering. The resulting paring symmetries 
in F are only ESE s-wave and OTE s-wave. It has been shown that OTE pairing se- 
riously influences on the density of states of quasiparticles. When OTE s-wave pairing 
dominates, quasiparticle density of state docs not have a gap-like structure around zero 
energy.®-'^'' In the extreme case, it has a zero energy peak (ZEP) structure. ^^'^^^ 

In §2.2 and §2.3, we restrict our discussion to ballistic transport regime. We will 
discuss relevant problems in diffusive regime in §2.4 and §2.5. 

2.2 N/S junction with ballistic system 

In this subsection, we consider superconducting proximity effect and resulting ubiq- 
uitous presence of odd-frequency pairing in N / S junctions. In general, physical quan- 
tities have spatial dependence near the interface. There arc two characteristic length 
scales: i) The inverse of the Fermi wave length and ii) coherence length ^ given by 
C — VfI where vp is the Fermi velocity and A is the energy gap (pair potential) in 
the superconductor. Here, we are interested in physical quantities varying slowly over 
the coherence length. In the present situation, it is appropriate to use quasiclassical 
Green's function method where high energy degrees of freedom are integrated out.^^"^^^ 
The quasiclassical Green's function is obtained from Gor'kov Green's function by inte- 
grating the energy and is concentrating on the low energy scale measured from the Fermi 
energy level. In this approximation, quasiparticles feel the pair potential depending on 
their direction of their motions on the Fermi surface. 

In the following, we consider two-dimensional N (a; < 0) / S (a; > 0) junctions in the 
ballistic limit (Fig. 2). For the bulk superconductor, we choose even- frequency ESE or 
ETO pairings. As regards the spin-triplet ETO pairing, we choose — Q for simplicity. 
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Fig. 2. (Color online) Schematic illustration of normal metal (semi-infinite) / superconductor (semi- 
infinitejjunction. A+,A_ denote pair potentials felt by right going (toward positive x) and left going 
(toward negative x) quasiparticles, respectively. 



There is no essential difference even if we consider spin-triplet pairing with 5*2; = 1. We 
assume a thin insulating barrier located at the N/S interface [x = 0) modeled by a 
delta functional potential H6{x). The reflection coefficient and transmission coefficient 
of the junction for the quasiparticle with the injection angle 6 is given hj R = Z'^/ (Z'^ + 
4cos^ &) and Tm = I — R, respectively, with Z = 2H/vf, where 9 (— 7r/2 < 9 < n/2) is 
measured from the normal to the interface. 



Here, we introduce quasiclassical retarded (advanced) Green's functions 
g^^^\e,9,r), and f2^^\£,9,r) where 9 is given by = A;ir(cos 6*, sin 6*). g^^^\e,9,r) 
{g^^^\e , 9 , r)) is the Green's function of quasiparticle at position r, energy e with 
injection angle 9 {tt — 9). f^^^\e,9,r) denotes anomalous Green's function which be- 
comes zero in the normal state. f^^^\e, 9, r) can be decomposed into f^^^\e, 9, r) and 
f2l^\e,9,r), 

fl^^\e, 9, r) = /fi^)(e, 9, r) - iffi^\e, 9, r) 

with 

/{i(£,0,r) = -At(-£,e,r), 

ff^{e,9,r) = if^{-e,9,r). (4) 

From this definition, it is clear that f^^\e^ 9, r) {f^^^\e, 9, r)) corresponds to the odd 
(even)-frequency pairing. Due to the translational invariance in the direction parallel to 
the interface x = 0, we can simply replace r by x. g^^^\e, 9,x) = g±}^\ f^^^\e, 9,x) = 
f^±^\ and f^^\e,9,x) = /^"^^ satisfy following so called the Eilenberger equation^^^ 



T^\v,.\^gl^^^ = [H,gl^^^] (5) 
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with 
and 

.R(A) _ R(A) . .R(A) . .R(A) . 

9± — 9± ^3 + /i± ri + /2± ^2- 

Here, (i = 1 — 3) is the Pauh matrix in the electron- hole space. If we write it explicitly 
decomposing each component, following three equations are derived 

I vfx I dxfi± = -2ief2± - 2A±{x)g±, (6) 
T« I vfx I dxf2± = 2i£/i±, (7) 

T« I ^^Fcc I d^g± = 2A±(x)/i±. (8) 

Here, we omit the superscripts R and A. All the Green's functions in eqs. (6)- (8) have the 
same analytical properties, i.e., R and A arc not mixed. In the above, vp^ = vp cos6 is 
the x-component of the Fermi velocity and A±{x) is the pair potential felt by ± direction 
in Fig. 2. In general, quasiparticle is scattered from electron(hole)-like quasiparticle to 
hole(electron)-like one by the pair potential. As seen from the Eilenberger equations, 
the scattering of quasiparticle by the pair potential is a driving force to produce pair 
amplitudes fi± and /2±. It should be remarked that A± and /2(i)± are different physical 
quantities. Even if A± is absent in a certain place x, /2(i)± may take a nonzero value 
by the proximity effect. A±{x) can be expressed by A±(x) = A{x)^±{9)<d{x).^^^ For s, 
dx2-y2, dxy, Px and Pj^- wave superconductors, $+(_)(^^) is given by 1 (1), cos 2^ (cos 2^), 
sin 2^ (—sin 26'), cos 6' (— cos^^), and sin^^ (sin^^), respectively. For sufficiently large 
magnitude of \x\ with x > 0, A(a;) converges to Aq, which is the bulk value of the pair 
potential. g±, fi± and /2± satisfy the normalization condition 

9l + fL + fL = 1- (9) 

We first note the general features of the pair amplitudes fi± and /2±. In normal state 
without A±, fi± = 0, /2± = 0, and g± = 1. This describes the quasiparticle without 
pair amplitudes. In bulk superconductors with A± — Aq for Aq > 0, we obtain fi± — 
since there is no spatial dependence of A±. For spin-singlet s-wave superconductor, 
g± = £/v/£^^^, /2± = Ao/^/Al - e'^. At £ = 0, ^± = and /2± = 1. This means 
the absence of quasiparticle and only the presence of Cooper pair at £ = 0. 

In order to discuss the parity with respect to the frequency much more clearly. 
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it is useful to use the Matsubara frequency representation: The real energy e is re- 
placed by iujn where the Matsubara frequency a;„ is given by a;„ = 27rT(n + 1/2) 
with temperature T and integer n. Here, we pay attention to the ioUn and 9 de- 
pendences of fi± = fi±{iun,0,x), f2± = f2±{ioJn,0,x). FoT X = oo, cqs. (6)-(8) 



lead to f2±{iuJn,9,oo) = Ao$±/a/A^$|+cj2, f2±{iuJn,0,x) = f2±{-iu}n,9,x) and 
fi±[iLJn,9,x) = 0. This means that only even-frequency pairing is possible in the 
bulk. Although fi±{iujn,9,x) is absent in the bulk, it has a nonzero value due to the 
spatial change of the pair potential A±. As seen from eqs. (6)-(8), fi±{iujn,9,x) — 
—fi±{—iujn,9,x) and f2±{iujn,9,x) — f2±{—iujn,9,x) are always satisfied. This means 
that the symmetry of ^) belongs to odd- frequency pairing. Furthermore, 

as seen from eqs. (6)- (8), 

v%dlfi± - 4(^^ + Al{x))h± ± 2z I I [d,A^{x)]g± = (10) 

is satisfied in the Matsubara representation. It is clear that the spatial dependence of 
the pair potential A.j-{x) is a source of the generation of odd frequency pair amplitude 
fi±- 

We discuss the parity of fi±{i(^n, 9, x) and f2±{ioOn, 9, x). ^±{—9) = ^z^{9) is satisfied 
for an even-parity superconductor, while ^±{—9) = —^^{9) is satisfied for an odd-parity 
one. It follows from eqs. (6)-(8) that fi±{iujn, 9, x) = —fi^iiun, —9, x) and f2±{i^m 9) = 
f-2^ii^m—9,x) for an even-parity superconductor and fi±{iujn,9,x) — f izf {iuJn, ~9,x) 
and f2±{iujn,9,x) — —f2^{i^n-i for an odd-parity superconductor. Note that the 

parity of the odd-frequency pair amplitude fi±{iujn,9,x) is different from that of the 
bulk superconductor. 

Next, we discuss the values of the pair amplitudes at the interface x — 0. For this 
purpose, it is convenient to express the above anomalous Green's function as^^^ 

_±iCF± + D±) _ D±-F^ 

l-D^F^ ' ^^^-1-L>±F±' ^^^^ 

where D± and F± satisfy the Eilenberger equations in the Riccati parameterization^^^ 

vfAD± - -A±(x)(l - Dl) + 2uJr,D± (12) 

VFxdA = -A±(a;)(l - F|) - 2a;„F±. (13) 

Since the interface is flat, F± — —RD^ for a;„ > and F± — —R~^D^ for Un < hold 
at x — 0^^) with reflection coefficient R. If we denote a;„ dependences of D± and F± 
explicitly, D±(a;„) = 1/ D±{—uJn) and F±{ujn) — 1/ F±{—uJn) are satisfled. 



J. Phys. Soc. Jpn. 



FULL PAPERS 



We concentrate on two extreme cases with (A) ^+{0) — ^-{0) and (B) ^+{0) — 
— $_(^). In the first case, there is no sign change of the pair potential felt by the 
quasiparticle reflected at the interface. Then, ABS is absent. On the other hand, in the 

second case, due to the sign change of the pair potential,^®-' ABS is generated near the 
interface. By using Riccati parameters defined above, it is easy to express that 



h±^±i{l-R)Dj{l + RDl), f2±^{l + R)Dj{l + RDl), a;„ > 

h± = ±z{R-l)D+/{R + Dl), f2± = {l + R)D+/{R + Dl), < (14) 



/i± = i(l + R)D+/{1 - RDl), h± = ±(1 - R)D+/{1 - RDl), a;„ > 

h±^i{l + R)D+/{R-Dl), f2±^^{l-R)D+/{R-Dl), a;„ < (15) 



for the case (B), respectively. In the low transparent limit with — > (i? — > 1), 
only f2± is nonzero for the case (A) and fi± is nonzero for the case (B) where Tm is a 
transmission coefficient. However, even in the case (A), the mixture of odd and even- 
frequency pairings occur at the interface for general R. The underlying physics behind 
these results can be qualitatively explained as follows. Due to the breakdown of the 
translational invariance, even- and odd-parity pairings are coupled near the interface. 
To be consistent with the Fermi-Dirac statistics, the interface-induced pair amplitude 
/i± should be odd in frequency with odd(even)-parity, when the pair potential has an 
even(odd)-parity.^°"^^) Note that /i± is pure imaginary and /2± is a real number, and 
there is a ±7r/2 phase between them. This twist of the phase is forced by the coexistence 
of the even- and odd-frequency components near the interface so as to be compatible 
with the time- reversal symmetry in the bulk superconductor. ^'^^ We summarize above 
results in Table II. 

Finally, we comment that analytical solution of Eilenberger equation is possible in a 
special case for fixed 9. If we choose spatial dependence of A+(x) as 



for the case (A) and 




A_(a;) = Ao tanh(a;/eo) 



(16) 



with a certain constant we obtain following results. 
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Table II. Paring symmetry in the bulk and that at the N/S interface. The pair amplitude in the 
bracket is suppressed in the low transparency limit. 



/2+ = -/2- = , ,2l.°A2 tanh(^). (17) 

It is evident that odd-frequency pair amplitude /i± is localized near a; = similar to 
zero energy quasiparticle state as seen from g±. Similar spatial dependence is realized 
in spin-singlet dxy-waye and spin-triplet Px-waye superconductor for low transparent 
limit. 

In §2.2.1 and §2.2.2, we illustrate these results by numerical calculations. In order 
to understand the angular momentum dependence of pair amplitudes, we define /i, /2 
as 

? J f^2)+ii^n,0,x), for -7r/2 < ^ < 7r/2, 

/l(2) = \ [^°) 

[ fi{2)-iicOn, 7r-e,x), for ti/2<9 < 37r/2. 
We decompose /i(2) into various angular momentum components. 

2.2.1 Proximity effect in spin-singlet s-wave superconductors 

We choose conventional spin-singlet s-wave superconductor as a typical example 
of ESE pairing. As seen from Table 11(1-1), both ESE and OSO pairing exist at the 
N/S interface. We plot in Fig. 3 the s-wave component of /2 Es{x) and the pj;-wave 
component of /i Op^{x). 

As seen from Fig. 3(a), the pair potential A.{x) is suppressed for fully transparent 
case (a) (T^ = 1, -R = 0) due to the penetration of Cooper pair into normal metal. At 
the same time, odd-frequency pair amplitude Op^{x) is enhanced at the interface. On 
the other hand, for low transparent case (b) (T^ < 1), Op^{x) is almost zero and A{x) 
remains constant up to the interface. 
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Fig. 3. (Color online) Spatial dependence of the pair potential for spin-singlet s-wave superconduc- 
tor. Distance x is normalized by ^ = vf/Aq. Pair potential normalized by its value in the bulk (solid 
line). Even frequency pair amplitude Ea{x) (dotted line), and odd frequency pair amplitude Op^{x) 
(dashed line). (a)Fully transparent case with Z = and (b)Low transparent case with Z = 5. a;„ = ttT 
and T — O.OSTp. Tq is the transition temperature of bulk superconductor. R = Z"^ /{Z"^ + 4cos^ 6). 



(a) (b) LDOS 



Normal metal Superconductor 




Fig. 4. (Color online) (a)Formation of standing wave by the closed loop of electron and hole. (b)Local 
density of state of quasiparticle p{e), 



In order to highlight the effect of the induced odd-frequency pairing more clearly, 
here we consider N/S junction with finite length of N^^^ [see Fig. 4 (a)]. We assume 
that the transparency at the interface is unity. In the present system, it is known 
that ABS with nonzero energy is formed''^-' as a standing wave which is generated by 
the interference between electron and hole in N. As seen in Fig. 4 (b), there are many 
peaks in local density of state (LDOS), which have been observed as Rowell-McMillan 
oscillation in tunneling spectroscopy. Here, we focus on this well-known phenomena 
in the light of odd-frequency pairing. Let us consider the pair amplitudes in the normal 
metal which are decomposed into components with different injection angle 9. We will 
define the following spatially averaged odd-frequency pair amplitude F[^\e,6) and 
even-frequency pair amplitude F2'^\e,6) 

Fi^,lie, 6) = ^ f,%^^{e, 9, x)dx. (19) 
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L=5Lo 




Fig. 5. Ratio of the spatially averaged odd-frequency pair amplitude to that of even-frequency one 
i?odd/even is plotted as a function oi e for = 0. Z = 0, L = 5Lo, and Lq = vf/{t^Tc). 



Then we introduce the ratio of F^^\e, 9) to F2l\e^ 9) as 



AN), 



R, 



odd/even 



tan 



/2eL 



(20) 



As seen in Fig. 5, in a certain energy regime, the magnitude of the odd-frequency 
pair amphtude exceeds over that of even-frequency one. It is remarkable that odd- 
frequency pairing can become dominant even in simple normal metal / spin-singlct 
s-wave superconductor junctions for high transparency at the interface. In the extreme 
case with L ^ Lq for fully transparent N/S interface, ABS levels with nonzero energy 
are given by™'^^-* 



2L 



'-{n + 1/2), 0,1,2, 



(21) 



It should be noted that just at e = £„, F2^\e,9) vanishes and -Rodd/cvcn diverges. 
Thus, the well-known quasiparticle bound states in N/S junction with conventional 
spin-singlet s-wave superconductor can be expressed by the generation of odd-frequency 
pairing. 

For low transparency limit, i.e., — > 0, the present ABS disappears since the 
couphng between N and S is weakened. Thus, the present OSO state is sensitive to the 
transparency at the interface. (We represent such a state in the braket in Table II.) 

In the next subsection, we will focus on spin-triplet p^^-wave superconductor, where 
odd-frequency pairing exceeds over even-frequency one at the interface independent of 
the transparency. 
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(b) low transparency 
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Fig. 6. (Color online) Spatial dependence of the pair potential for spin-triplet px-wave superconduc- 
tor similar to Fig. 3. Distance x is normalized by ^ = vp/Aq. Pair potential normalized by its value 
in the bulk (solid line), Even frequency pair amplitude Ep^{x) (dotted line), and odd frequency pair 
amplitude Os{x) (dashed line). (a)Fully transparent case with Z — (b)Low transparent case with 
Z = b. w„ = ttT and T = 0.05Tc. 



2.2.2 Andreev bound state in unconventional superconductor 

In this subsection, we show that the mid gap ZEABS specific to unconventional 
(anisotropic) superconductor is expressed by odd-frequency pairing. We focus on the 
case with $_|_$_ < where mid gap ABS is generated. In this discussed in 

§2.2, eq.(15) is satisfied by using Riccati parameters. In the low transparency limit 
with Tr„ — 7- (i? — 7- 1), only the fi± is nonzero. Namely, only the odd-frequency pair 
amplitude exists at the interface for — !■ 0. First, we focus on the spatial dependence 
of the pair potential and pair amplitudes for spin-triplet p^-wave superconductor. We 
plot in Fig. 6 the p^^-wave component of ETO pairing Ep^ (x) and the s-wave component 
of OTE pairing Os{x). 

For fully transparent junctions, similar to spin-singlet s-wave junction case (Fig. 
3(a)), odd-frequency component of the pair amplitude Os{x) becomes of the same or- 
der of even-frequency one Ep^{x) as shown in Fig. 6(a). For low transparent junction, 
Os{x) largely exceeds over Ep^{x) near the interface while A(x) is suppressed [Fig. 
6(b)]. The corresponding LDOS p{e) at x = has a ZEP for spin-triplet p^^-wave su- 
perconductor junction (curve (a) in Fig. 7), while p{e) has a gapped line shape for 
spin-singlet s-wave superconductor junction. Similar ZEP also appears for spin-singlet 
dxy-'wave superconductor. ZEP of p{e) ubiquitously appears for unconventional super- 
conductor junctions, where the pair potential changes sign on the Fermi surface.^^'^^'^^^ 
Experimentally, ZBCP by this ZEP has been observed in many unconventional super- 
conductors including cuprate.^^'^^''^^"^^'^^"^^-' It is noted that in the presence of the ZEP 
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Fig. 7. (Color online) Quasiparticle density of state p{e) normalized by its value in the normal state 
at the N/S interface for Z = 5. (a)Spin-triplet Pa;-wave superconductor junction, (b) Spin-singlet s-wave 
superconductor junction. 



originating from ABS, the magnitude of the odd-frequency pairing is enhanced. 

Hereafter, we discuss more details about the odd-frequency pairing in the low trans- 
parency limit. As shown in the last paragraph, only fi± exists for i? — )■ 1 with < 0. 
In this limit, /2± becomes zero and fi± is given by fi± = 2iD^/{l — -D^). Furthermore, 
ABS is expressed only by the induced odd-frequency pairing. Then, /i in eq.(18) is 
given by 

^ r ^l£^ -n/2<9<n/2 

\ -2iD-{-K-e) _/0//)/Q_/o ^ ' 

{ 1-D-il-ey vr/2<0<37r/2 
with D+{e) = -D^{e) and D^{7t - 6) = D+{0). Here, we explicitly write the 6 de- 
pendence of D^. When the parity of bulk superconductor is even [odd], D^{9) = 
D+{tt + 6)= D_{-e) [D+{9) = -D+irr + 6) = is satisfied. Thus, A has an 

odd-parity (even-parity), when the bulk paring symmetry is even-parity (odd-parity). 
The difference in the parity of the induced odd-frequency pairing results in a serious 
difference when we consider proximity effects into DN attached to superconductor.^*^'™) 
In DN, only s-wave even parity pairing is possible. Thus, pair amplitudes with angular 
momentum larger than cannot penetrate into DN. This implies that odd-frequency 
pairing generated from dxy-wave superconductor cannot enter into DN^*^'*^^) since it is 
expressed by OSO pairing. On the other hand, for ]9^.-wave superconductor, ABS can 
enter into DN®^"*^^) since it is expressed by OTE including s-wave channel. 

If we neglect the spatial change of the pair potential near the interface, it is possible 
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to obtain and / analytically. Since D+ is give by A+ / {ujn + A^) , /i is given 

by 



»A+(e) 



-7r/2 < ^ < 7r/2 



For spin-triplet p^-wave bulk superconductor, the resulting /i is 

h = cost' , 

while for spin-singlet dxy-wave bulk superconductor, the resulting /i is 

/i = - — - I sin 26* I sgn(sin^). 

Reflecting on the ZEABS, /i is proportional to the inverse of Un for both cases. The 
difference of the parity of /i between two cases appears in the topological features of 
wave function of two ABSs as will be discussed in 53.3. 



2.3 Odd-frequency pairing in vortex core 

In this subsection, we discuss the pairing symmetry in vortex core.^^^ The study of 
the mixed state in type-II superconductors, where magnetic flux enters a sample in the 
form of quantized vortices, has a long history. It is known that the ABS is generated 
in the Abrikosov vortex core due to the spatial structure of the pair potential. ^^"^^^ 
One of the manifestations of the bound states is the enhancement of LDOS in the 
core, observable as a ZBCP by scanning tunnehng spectroscopy (STM).^°'^^) Since an 
Abrikosov vortex breaks translational symmetry in a superconductor, we can expect the 
emergence of an odd-frequency pairing state around the vortex even in a conventional 
spin-singlet s-wave superconductor. 

Here, based on the quasiclassical theory, we clarify a pairing symmetry in the vortex 
core. The electronic structure at the core of a single Abrikosov vortex in a ballistic super- 
conductor is described by the quasiclassical Eilenberger equations.^®'^^'^^^ We assume 
that the pairing symmetry of bulk superconductor is ESE or ETO. Along a trajectory 
r{x') — ro + x' Vp with unit vector Vp parallel to the Fermi velocity Vp, the Eilenberger 
equations are generally represented in a 4x4 matrix.^^^ For a singlet (triplet) supercon- 
ductor with Sz — 0, these equations are reduced to the set of two decoupled differential 
equations of the Riccati type for the functions a{x') and b{x')^^'^^^ 

VFd^,a{x') [2un + A*a(x')] a{x') - A 0, 
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vpd^'bix') - [2ujn + Ab{x')] b{x') + A* = (24) 

where cOn is the Matsubara frequency. Quasiparticle Green's function g and anomalous 
Green's function / are expressed by a and b as g = {l — ab)/{l + ab) and f = 2a/ {1 + ab), 
respectively. For simple case of a cylindrical Fermi surface, the Fermi velocity can be 
written as Vp — vf{cos9, sin^). We choose the following form of the pair potential: 

A = A(r, 9) = Ao$(6')F(r) Gxp{imip) 

with r = ^/x^~+y^, exp{i(p) — {x + iy)/r and integer m. The center of a vortex is 
situated at x — y — and exp{im(p) is the phase factor which originates from the 
vortex. We consider axially symmetric cores. For the calculation of the normalized local 
DOS by its value in the normal state, the quasiclassical propagator has to be integrated 
over the angle 9 which defines the direction of the Fermi velocity. The normalized LDOS 
in terms of functions a and b is given by 

P(ro,s)^ / -Re — — , (25) 

Jo 27r ll + ab\^^^^^^s 

where e denotes the quasiparticle energy with respect to the Fermi level and d is an 
effective broadening parameter of quasiparticle energy level. In the actual numerical 
calculations, we will fix this value as 6 = 0.1 Aq. 

First, we show the general property of the symmetry at the vortex center inde- 
pendent of the detailed spatial dependence of A. Let us consider a trajectory passing 
through the center of the vortex. By setting x' — at the center of the vortex, we ob- 
tain b{x',iujn) — —l/a{—x', —iujn) from the Eilenberger equations with odd integer m. 
Similarly, we obtain b{x', ioun) — l/a{—x', —iuin) for even integer m. Thus, at the vortex 
center x' — 0, we get fiioJn) = ~f{~i'^n) in the former case, while fiioJn) = f{—i'^n) 
in the latter. Since we do not consider the Zeeman effect, spin flip does not occur in 
the present system. Thus, we can summarize pairing symmetry at the vortex center in 
Table III based on the broken inversion symmetry and Fermi-Dirac statistics. 

For conventional ESE s-wave case, there have been several studies of multi-flux state 
with m > 1 100' 101) in the context of superconducting quantum dot where normal nano- 
scale region is surrounded by superconductor. It has been shown that ZEP in the 
DOS only appears for odd number m at the vortex center. ^°°'^°^) This statement is 
consistent with our present results for the conventional spin-singlet s-wave case, where 
odd-frequency pairing is generated only for odd integer m. 

In general, the most stable vorticity is | m |= 1 realized in Abrikosov vortex. As 
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Table III. Pahing symmetry in the vortex state. 




Fig. 8. (Color online) (a) Normalized LDOS around the vortex at e = 0. The center of the vortex 
is situated at x — y ~ 0. (b) Spatial dependencies of various pairing components at e = 0. Spatial 
dependencies of (b) ESE component and (c) OSO component. [Reproduced from Fig. 1 of Phys. Rev. 
B 78 012508 (2008) by T. Yokoyama, et. al.] 

a symmetry of the bulk superconductor, we choose ESE s-wave pairing. Also, spatial 
dependence of the gap is chosen as F(r) = tanh(r/,^), where ^ = vf/Aq is the coher- 
ence length. Due to the broken translational symmetry of the system, various pairing 
symmetries are mixed near the vortex within the spin-singlet pairings. We decompose 
anomalous Green's function f^{e, 9) into various angular momentum components at 
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the center of the core, and f^{e, 9) can be given by 

/«(£,^)= fi.^M^hO). (26) 

;z=o,±i,±2,... 

All the above pairing components fi^ are spin singlet. Spatial dependence of the LDOS 
around the vortex at £ = is shown in Fig. 8 (a). It has a peak just at the core center. 
The spatial dependencies of decomposed anomalous Green's function f]^ at £ = are 
shown in Figs. 8 (b) and (c). Interestingly, only OSO pairing component Re(/i) is 
nonzero at the center of the vortex core. Thus, we see that anomalous Green's function 
at the core center has chiral p-wave symmetry. 

It is interesting to consider the bulk chiral superconductor with the an- 
gular momentum I. The corresponding pair potential is given by A = 
Ao$(^)F(r) exp(i/2.^) exp(im(^). At the center of the core, it is shown analytically that 
the angular momentum of the pair amplitude is + m. Here, we choose Abrikosov 
vortex in bulk chiral p-wave superconductor with = 1. As seen from Table 111(3), 
symmetry at the center of core is OTE. The angular momentum of the pair amplitude 
depends on whether chirality and vorticity is antiparallel or parallel. For the antiparallel 
case {m — —1), the OTE s-wave pairing is generated at the vortex core,^°^^ while for the 
parallel case (m = 1), the OTE d-wave symmetry is realized at the core. Since s-wave 
pair amplitude is robust against impurity scattering,^°) it is useful to consider impurity 
scattering effect on the vortex core state (odd-frequency pairing) to distinguish two 
states. We introduce this effect within the Born approximation, where we denote the 
impurity scattering rate in the normal state F with F = 1/ (2r) with mean free path 
VpT. The resulting LDOS of anti-parallel and parallel vortex case is shown in Fig. 9. 
The ZEP is robust with the increase of F for anti-parallel vortex case while it is fragile 
against F for parallel vortex case.^^^"^"^^ 

The robustness of the ZEP at the vortex center against nonmagnetic impurities for 
antiparallel case originates from the odd-frequency s-wave pair amplitude. In the actual 
sample of spin-triplet chiral p- wave superconductor, degenerated chiral state, i.e., 1^ — 1 
and Iz — —1, form a domain structure. Near the antiparallel vortex core, strong ZEP of 
LDOS is expected, while near the parallel vortex core, weak ZEP appears. Accordingly, 
measurements of the ZEP of LDOS in the presence of impurities reflects the detection 
of the symmetry of odd-frequency pair amphtude at the center of the core as shown in 
Fig. 10. 

Extension of the calculation in the Abrikosov vortex lattice has been done by 
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(a) Antiparallel vortex (b) Parallel vortex 
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Fig. 9. (Color online) The LDOS at the vortex center for the antiparallel (a) and (b) the parallel 
vortex. The dashed and solid lines are plots for F = O.IAq and O.SAq, respectively. [Reproduced from 
Fig. 3 of Phys. Rev. Lett. 102 117003(2009) by Y. Tanuma et. al. ] 



(a) Antiparallel chiral p-wave votex (b) Parallel chiral p-wave votex 




Fig. 10. (Color online) Schematic illustration of odd-w pair amplitude for (a) the antiparallel and 
(b) the parallel vortex in chiral p-wave SCs. The arrows represent the phase rotation. [Reproduced 
from Fig. 4 of Phys. Rev. Lett. 102 117003(2009) by Y. Tanuma et. al. ] 

Yokoyama et. al.^^'^^ They have found that only odd- frequency pairing exists at core 
centers. 

2.4 Anomalous proximity effect in spin-triplet superconductor junctions 

In this subsection, we discuss the anomalous proximity effect specific to spin-triplet 
superconductor junctions. We consider diffusive normal metal(DN) (0 < x < L) / 
superconductor (S) (x > L) junction where the length of DN is L as shown in Fig. 11. 
We assume that the DN is attached to the normal electrode at x = 0. The interface 
between the DN and the superconductor (S) at x = L has a resistance Rb and the 
N/DN interface at x = has a resistance Rb>. We also denote resistance in DN as Rd- 
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Fig. 11. Schematic illustration of the model of normal electrode attached to diffusive normal metal 
(DN) / superconductor (semi-infinitejjunction. [Revised from Fig.l of Phys. Rev. B 72 140503(R) 
(2005) by Y. Tanaka et. al.] 



For Rf)/ = oo, the present model is reduced to the DN/S bilayer with vacuum at the 
DN free surface. 

We assume that pairing symmetry in S is ESE or ETO in the balhstic regime. We 
choose Sz = for ETO case. In the DN region, only the s-wave pairing which is robust 
against impurity can survive in the limit of strong disorder. Since there is no spin-flip 
in the present junctions, spin structure of pair amplitude in DN is the same as that in 
bulk. Then, the pairing symmetry in the DN can be derived in accordance with Fermi- 
Dirac statistics. If the pairing symmetry in the bulk is ESE, the resulting symmetry in 
DN is ESE s-wave. On the other hand, if the pairing symmetry in the bulk is ETO, 
the resulting symmetry in DN is OTE s-wave. In DN, it is convenient to use Usadel 
equation which is derived from Eilenberger equation in the diffusive limit after angular 
average of the direction of motion. The Green's functions in Usadel equation do not 
have a 6 dependence anymore due to the isotropization by impurity scattering. Here, 
in order to calculate LDOS in the DN region, we use retarded Green's function. Due 
to the fermion's commutation relation and complex conjugation of Green's function, 
anomalous Green's function should satisfy f^2){£,x) = [/i(2)(^' ^)]* where the suffix 1 
and 2 denote odd and even-frequency component. f(^{e,x) = —[ff'{—e,x)]* is satisfied 
by the relation f^{e,x) = —f^{—e,x). On the other hand, f^{e,x) = [f,^{—e,x)]* is 
satisfied by the relation f^{e,x) = f,^{—e,x). Usadel equation^^^^^ in DN is given by 

D-^ + 2iesmC = 0, (27) 

with diffusion constant D. g^{6,x) and f.^{e,x) are parameterized by g^{e,x) = cos(, 
f^{e, x) = sin (, and f^{e, x) = 0, when the bulk superconductor has an ESE symmetry. 
On the other hand, g^{e,x) = cos(, fi{e,x) = sin(, and = are satisfied 
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when the bulk superconductor has an ETO symmetry. The boundary condition of ( is 
given by^°"^^) 



x=L 



(28) 



p ^ 2Ti{fs COS Cl - gs sing) ^29) 
1 2-Ti+Ti(cosCL5sH-sina/5) ^ ' 



at DN/S interface {x = L) and 



(^^) = „ f'Z'" . . (30) 



^=0 ^b' 2-T2+T2 cos Co ' 

at normal electrode /DN interface {x = 0). Here, (l = C \x=l, and Co = C U=o- (• • •) 
denotes the angular average of the injection angle 9 of quasiparticles, 

fTr/2 



JZ'^,d0 COS 0F^2){e) 



i^mio)) = _ "V , (31) 



_ 4cos^6' _ 4cos^6' 

' ^ Z2 + 4cos2^' Z'2 + 4cos2^' ^ ^ 

Ti^2 denotes the transmission coefficient at the interface with barrier parameters Z 
and Z'. gs and fs are determined by the pair potential in bulk superconductor with 
gs = {g+ + g-)/{l + g+g- + /+/_), g± = e/V^^^, /± = A±/yAl^, and 
A± = Ao$±. For even-parity bulk superconductor, fs — (/+ + /-)/(! + 9+9- + /+/-) 
and for odd-parity superconductor fs — i{f+g-—f-g+)/{l+g+g-+f+f-).^'^'^^'' Here, we 
neglect the spatial dependence of A±. By solving, above Usadel equation with proper 
boundary condition we can determine the pair amplitude in DN and the resulting 
density of state. 
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Table IV. Paring symmetry in the bulk, at the interface, and in DN. 
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The obtained results are summarized in Table IV. LDOS p{e) normalized by its value 
in the normal state is plotted in Fig. 12. The pair amphtude /i and /2 at the interface 
of S side discussed in §2.2 can not enter into DN if their angular momentum I is not 
zero. When the symmetry of the bulk superconductor is ESE s-wave (1-1), fl^{s, 6) = 
and f-^{e, 9) are satisfied. This means that ESE s-wave pair is induced in DN. The 
OSO pair amplitude induced at the interface can not enter into DN due to the odd- 
parity. The LDOS has a gap like structure as shown in Fig. 12(a). This is a standard 
proximity effect known for the mesoscopic superconducting systems. The order 
of the magnitude of the energy gap is Thouless energy Eth — D/L"^ in DN. When the 
symmetry of the bulk superconductor is ESE dxy-wme (1-2), the resulting ( is zero 
and f^{e, 9) — f^{e, x) — inside DN. Then the resulting p{e) is unity as shown in 
Fig. 12(b). The OSO pair amplitude at the interface can not enter into DN. Also, the 
subdominant ESE state can not enter since its angular momentum / is / > 2. This 
is the reason for the absence of the proximity effect into DN. In other words, mid 
gap ABS generated at the interface can not penetrate into DN for bulk spin-singlet 
(ij^y-wave superconductor. For ETO state with py-waye pairing (2-1), ( is zero and 
f^{e, 9) — f^{e, 9) — inside DN. Then the resulting p{e) is unity as shown in Fig. 
12(c). In this case, pair amplitudes at the S side of the interface have angular momenta 
with I > 0. Proximity effect specific to spin-triplet superconductor appears for ETO 
Px-wme bulk superconductor case (2-2). The resulting C and f^{e,x) are nonzero m 
DN with f'^{e, x) = 0. It is remarkable that p{e) has a sharp ZEP. This is a new type of 
proximity effect. The present anomalous proximity effect is generated by odd- frequency 
pairing. '^^^ The mid gap ABS at the interface can penetrate into DN.*^^'*^"^) In this case, 
mid gap ABS can be interpreted as an OTE pairing which has an s-wave component.™^ 
Then, the proximity effect into DN becomes possible. 

Although we have shown anomalous proximity effect for ETO p^-wave superconduc- 
tor, it appears for generic ETO superconductor as far as the the interface induced OTE 
pair amplitude has an s-wave component. ^^'^^^ ETO p^-wave case is exceptional since 
OTE s-wave component is canceled by the angular integral of 9 in the boundary con- 
dition. One of the strong candidate material of ETO p-wave superconductor is realized 
in chiral p-wave superconductor Sr2Ru04 where ZBCP originating from ABS^^'^^^^ has 
been observed. ^^"^^^ In this case, ZEP of p{e) in DN by the anomalous proximity effect 
appears. In the actual calculation, we must extend boundary condition at DN/S inter- 
face.^°^ In Fig. 13, we plot LDOS p{e) for various parameters. The height of the ZEP 
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Fig. 12. LDOS p{e) of quasiparticle inside the diffusive normal metal at x = L/2. Symmetry in the 
bulk superconductor is (a)ESE s-wave, (b)ESE rfxj/-wave, (c)ETO Py-wave, and (d)ETO Px-wacve. 




Fig. 13. Similar plot to Fig. 12 of p{e) inside the diffusive normal metal at a; = L/2 for ETO chiral 
p-wave superconductor. Z ^ Z' = 1, Rd/Rb = 1 and Erh = 0.25Ao. {a)Rd/Rv = 0.01, {h)Rd/Rb' = 1, 
and {c)Rd/Rb' = 100. [Revised from Fig.2(a) of Phys. Rev. B 72 140503(R) (2005) by Y. Tanaka et. 
al] 



of LDOS increases with the decrease of Rd/Rb' since the odd pair amphtude fi{s, x) is 
more strongly confined within DN for large magnitude of Ry. 

In order to detect anomalous proximity effect by odd-frequency pairing, we propose 
an experimental setup using scanning tunneling spectroscopy (STS) as shown in Fig. 14. 
The ZEP of LDOS can be detected by scanning tunneling spectroscopy (STS) and we 
can distinguish anomalous proximity effect by OTE pairing from conventional proximity 
effect by ESE pairing. However, to fabricate well-controlled surface of junctions available 
for STS is not so easy at present. In order to detect anomalous proximity effect, Asano et. 
al proposed a T-shaed junction®^) which is more accessible for the fabrication as shown 
in Fig. 15. In this system, the proximity effect from the superconductor modifies the 
conductance between two electrodes depending remarkably on the pairing symmetry: 
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(a) Anomalous proximity (b) Conventional proximity 

STS STS 
▼ ▼ 

^ OTE ^ ETO y y ESE ^ ESE ^ 
DN Superconductor DN Superconductor 

Fig. 14. Schematic illustration of an experimental setup by scanning tunneling microscope (STM) in 
diffusive normal metal (DN) attached to a superconductor. (a)Anomalous proximity effect in DN/ETO 
superconductor junction, (b) Conventional proximity effect in DN/ETO superconductor junction. The 
pair amplitude has an OTE and ESE symmetry for anomalous case (a) and conventional case (b), 
respectively. Only for anomalous proximity case (a), the LDOS has a zero energy peak. 



ESE ESE ETC 




Fig. 15. (Color online) Schematic illustration of an experimental setup by T-shaped proximity struc- 
ture. We consider spin-singlet s-wave, d-wave and spin-triplet p-wave pair superconductors. [Repro- 
duced from Fig.l of Phys. Rev. Lett. 99 067005 (2007) by Y. Asano et. al\ 



spin-singlet or spin-triplet. Only for spin-triplet pairing, OTE state can be generated 
and conductance has a zero bias peak. 

The anomalous proximity effect induces an peculiar response of DN to external 
magnetic field. When a magnetic field is applied parallel to the interface, the magnetic 
field in the DN behaves as if(x) ~ exp(— x/A(x)) with the local penetration depth A(x), 
which is given by 

1 _ TY..jl{^^n.x) 

using Matsubara representation as shown in eq. (18) of ref.^^^-* or eq. (2.5) of ref.^^^-* 
with Ao ^ = 32n^e^N{0)D. In Fig. 16, the averaged value of A^, (A^^ = L/ ^) is 
plotted as a function of temperatures (T), where Tc is the transition temperature of bulk 
superconductor. As is shown in Fig. 16 by curves c and d, A^^ > for corresponding ESE 
superconductor junctions. Thus Xav is a real number and a magnetic field is screened by 
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the usual Meissner effect in tfie DN irrespective of the fact whether bulk superconductor 
keeps T symmetry or not. On the other hand, in ETO superconductor junctions, we 
find A^^ < as shown by curves a and b. Therefore, the A^^ becomes a purely imaginary 
number in the anomalous proximity case. This is the consequence of the fact that the 
pair amplitude fi{iuJn,x) is purely imaginary. The negative value of A^^ means the 
generation of negative superfiuid density locally. It is a novel feature of the anomalous 
proximity effect that the applied magnetic field is not screened in the DN region and 
is screened only by the S region. The magnetic field can spatially oscillate in DN and 
paramagnetic state becomes possible. The paramagnetic Meissner effect in the surface 
of high Tc cuprate has been studied^^^^ and anomalous temperature dependence of 
the penetration depth has been reported. ^^^'^^^^ Recently, Asano et. al have calculated 
dynamical response of OTE pairing to the electromagnetic field. ^^^^ It has been found 
that the surface impedance {Z = R — iX) of a DN has anomalous features where R 
and Z are resistance and reactance, respectively. In contrast to the standard relation 
[R « X), an anomalous relation i? > X is satisfied at low temperatures. 

Before closing this subsection, it is noted that purely odd-frequency pairing is real- 
ized in DN for DN/ETO superconductor junctions. The present odd- frequency pairing 
in DN arouses plenty of anomalous quantum phenomena including enhanced Joseph- 
son current in S/DN/S junctions. Since there are many new quantum phenomena 
relevant to anomalous proximity effect via odd-frequency pairing, we hope it will be 
verified by experiments in Sr2Ru04 junctions. ^'^'^^'^^^^ 



0.02 




T/Tc 1 



Fig. 16. Averaged value of the Meissner screening length A^^ is plotted for various superconductor 
with Z = 1, Rd/Rb = I, Z' = 1, Rd/Rb' = 1 and Eth = 0.25Ao, respectively. (a)ETO chiral p-wave, 
(b)ETO px -wave, (c)ESE s-wave and (d)ESE chiral d-wave superconductors. [Revised from Fig. 3 of 
Phys. Rev. B 72 140503(R) (2005) by Y. Tanaka et. al.] 
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2. 5 odd-frequency pairing in ferromagnet / superconductor junctions 

In this subsection, we discuss proximity effect in diffusive ferromagnet (DF)/ su- 
perconductor (S) junctions. The presence of odd-frequency pairing in DF/S junction 
has been originally proposed by Bergeret Volkov and Efetov.^^'^^^-' In their work, since 
inhomogeneous ferromagnet is assumed, equal-spin-triplet odd-frequency pairing is gen- 
erated by the spin flip of an electron. ^^^^ Then, so called long range proximity effect 
appears, where the characteristic length of proximity effect is not y^D/h but 
\J~DJT similar to the conventional proximity effect without exchange energy /i, where 
DF becomes diffusive non-magnetic metal denoted as DN in previous subsections. It is 
noted that there have been several relevant theoreticaP^^"^^^^ and experimentaP^^^^^^ 
works in DF/S junctions up to now. Although original idea by Bergeret et. al has as- 
sumed inhomogeneous magnetization, it is possible to generate odd-frequency pairing 
even if we assume uniform ferromagnet with constant h in DF. Here, we discuss the 
symmetry of pair amplitudes and resulting LDOS.^^^ We assume ESE superconduc- 
tor and ETO one with Sz = 0. In this case, straightforward extension of the previous 
subsection becomes possible. 

Before calculation of LDOS by Usadel equation, we discuss general properties. In 
DF, only even parity s-wave pairing symmetry is possible due to the impurity scattering. 
When the symmetry of bulk superconductor is ESE, original pairing symmetry in DF 
without h is ESE as discussed in previous subsection. By the nonzero h, OTE pairing 
is induced by the explicit spin rotational symmetry breaking. On the other hand, when 
the bulk symmetry is ETO, only OTE pairing is possible in DF without h. Then, ESE 
pairing is induced by h. Thus for generic case, both ESE and OTE pair amplitudes 
exist in DF. Hereafter, we assume that the magnitude of h is much smaller than Fermi 
energy in S and DF. In this case, it is possible to apply Usadel equation. The boundary 
condition of the Green's function discussed in previous section is also available now. In 
the present case, we must pay attention to the direction of Sz in the Green's function 
which we have not explicitly written in the previous sections. As in the case of §2.4, 
we neglect the spatial dependence of pair potential in the bulk. Then the bulk Green's 
function g± and f± are given as follows 

g+ = gf^^{e,e)^gf^^{e,e) 
g. = gf^^ie, 7^-9)= g^^ie, tt - 6) (34) 
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U^ff^^{e,9) f.^ff^^{s,n-9). (35) 

For ESE superconductor, f^^{e,9) = -f^^{e,9) is satisfied while f^^{e,9) = f^^{e,9) 
is satisfied for ETO one. g± and f± are given by g± — sj yje^ — and /± = 



A±/ — respectively, with A± = Ao$± as in the last subsection. In DF, C, 
in eq. (27) follows 

L'^ + 2i(£ + /i)sinC = 0. (36) 
with sin C, = f^ii^, x). On the other hand, sin ( with sin ( = f^^^is, x) satisfies 

D-4 + 2i(£-/i)sinC = (37) 

The boundary condition of C is given by eqs. (29) (30).^'''^^^ On the other hand, the 
boundary condition of ( is given by 



at the DF/S interface. Fi is given by 



x=L 



f^, (38) 



^ 2^l(/g^QSCL_-^5 SinC^) ^gg^ 

2 - Ti + ri(cos (l9s + sin Cl/s) 

for ETO superconductor and 

^ 2Ti(-/g cos Cl - sin (l) ^^^^ 
2 - Ti + Ti (cos Cigs - sin Cl/s) 
for ESE superconductor. At the interface between normal electrode and DF, 



{F2) ^ 2T2sinCo 



F2 = ^ z;' Z J - (41) 



.=0 ^6 2-T2 + T2COSC0 

Here, C,l — C, \x=l^ Co = C U=o- Even-frequency and odd-frequency pair amplitudes are 
given by 

/2^(£,a;) = (sinC-sinC)/2 (42) 

and 

//^(£,x) = (sinC + sinC)/2, (43) 

respectively. After solving Usadel equation, we can determine the symmetry of pair 
amplitudes in DF as shown in Table V. Proximity effect is absent for cases (1-2) and (2- 
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wave) 


No 


No 


(2-1) 


@ETO {py- 
wave) 


No 


No 


(2-2) 


ETO (p^- 
wave) 


OTE 


OTE + ESE 



Table V. Symmetries of Cooper pair in Ferromagnet (diffusive) / superconductor junctions. Tlie 
bold letter expresses tlic symmetry induced by the exchange field h in DF. 

1). For (1-1) and (2-2) cases, pairing symmetries in DF are ESE and OTE, respectively. 

In the following, we calculate f^{e, x), ff{e, x), and LDOS p{e). p{e) in DF is given by 

p(e) = ^ (Re cos C + Re cos C) (44) 

We consider a model similar to Fig. 11, where DN is replaced by DF. We assume ESE 
s-wave superconductor as a bulk state and choose Z — 3, Z' — 3, Eth = -D/L^ = O.lA, 
Rd/R'b — 0.1. Pair amplitudes at x — 0, i.e., the interface between normal electrode 
and DF, are plotted in Fig. 17. 

From the definition of Green's function, imaginary part of f^{£, x) and real part of 
ff{e,x) become zero at £ = for any h. For h = 0, f^{e,x) is zero and only f^{e,x) 
has a nonzero \ increase of h, the magnitude of /2 = f^{e,x) 

around e = is suppressed. On the other hand, the magnitude of the imaginary part of 
/i = f(^[e,x) is enhanced by h. The corresponding p{e) is shown in Fig. 18. For h = 0, 
p{e) has a mini gap in consistent with previous theory .^^^^ p(e) is sensitive to h. For 
h/A — 0.05, p{e) has a ZER The presence of ZEP of p{e) in F/S junction has been 
discussed in previous works.^2^'^^^-^^^) It has been clarified that for Eth ~ 2hRb/Rd, p{e) 
has a ZEP.^^^'^^^) This condition is consistent with the present choice of the parameters. 
As seen from Figs. 17 and 18, p{e — 0) is enhanced when the imaginary part of //^(s, x) 
has a large value (Fig. 17(d)). It is noted that odd-frequency pair amphtude plays a 
pivotal role for the generation of ZEP.^^^'^^^^ 

In Fig. 19, we study the crossover between ESE and OTE pairings. We plot Re/2 
and Im/i as a function of /i for e = at (a) x = 0, (b) a; = L/2 and (c) x = L in 
Fig. 19. Re/i increases from zero with h and it has a maximum at a certain value of h. 
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Fig. 17. (Color online) Real and Imaginary parts of pair amplitudes /2 = /2^(e, x) and /i = fi '{e, x) 
at a; = are plotted for DF/S junctions with Rd/Rb = 1 for various exchange energy h. The pairing 
symmetry of /2 and /i are ESE and OTE, respectively. [Reproduced from Fig. 2 of Phys. Rev. B, 
75,134510,(2007) by Yokoyama et. al] 




Fig. 18. (Color online) Normalized quasiparticle density of state pie) at x = for DF/S junctions is 
plotted as a function of e for Rd/Rb = 1- [Reproduced from Fig. 3 of Phys. Rev. B, 75,134510,(2007) 
by Yokoyama et al.] 



As shown in Fig. 19(a) and Fig. 19(b), if the value of h is larger than this value, OTE 
pairing becomes dominant. If we use sin(^(£:) = — sin(^*(— e) and cosC,{e) = cos(*{—e) 
with ( = ({e) and ( = C{e), following relations are satisfied in pair amplitudes in DF, 

/i = A^(£,x) = [sinC(e) + sinC(-£)]/2, (45) 

f2 = f2^{e,x) = [sinC(e)-sinC(-5)]/2. (46) 
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The ratio of /2 to /i at e = is given by 

/2 tanReC(O) 



(47) 

/i itanhlmC(O) ^ ' 

If |ReC(0)| < |ImC(0)| is satisfied, the crossover occurs and the OTE pair amphtude 
becomes dominant. The threshold h, where the present crossover occurs, is given by 



/;f2 
/ '. 


(a): 










(A) 








Xfl 


ic) - 





\> Q.I ^2 0.3 0.4 0.5 



Fig. 19. (Color online) Pair amplitudes with ESE symmetry /2 = /^ (e, a;) and that with OTE 
symmetry /i = f[^'{e,x) at e = are plotted as a function of h. (a) x = (normal electrode /DF 
interface), (b) x — L/2 (middle of DF). (c) x — L (DF/S interface). [Reproduced from Fig. 5 of Phys. 
Rev. B, 75,134510,(2007) by Yokoyama et. aL]j 



The generation of OTE state is possible without using DF region, i^^^^^^) hinder et 
al. has clarified that even in DN/S junctions, if the interface is spin active, OTE pairing 
can be induced in DN where the bulk symmetry is ESE s-wave.^^^'i^i-* In this case, also 
ZEP of LDOS can be expected. Recently, Yokoyama et. al has predicted the presence 
of anomalous Meissner effect is generated by the spin active interface. The magnetic 
susceptibility has a nonmonotonic temperature dependence accompanied by its sign 
change. Correspondingly, magnetic field and current density can spatially oscillate in 
DN.143) hope these features will be observed by experiments with yiiSR or microwave 
resonance. 

Proximity effect in fully polarized ferromagnet, i.e., half metal(HM), has become a 
hot topic now. Keizer et. al reported the existence of Josephson coupling in supercon- 
ductor/half metal/superconductor (S/HM/S) junctions, ^^^^ where symmetry of super- 
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conductor is conventional ESE s-wave. One of the possible pairing carrying a Josephson 
current in HM is a ETO p-wave pairing generated by the spin flip scattering at the in- 
terface. In real S/HM/S junctions, however, half metals are close to the dirty limit in 
the diffusive transport regime. Then more promising symmetry which can carry Joseph- 
son current is OTE s-wave pairing. Stimulated by the experiment by Keizer et. al, there 
have been several theoretical proposals^^^"^^^-* and an experimental report. Since the 
magnitude of the exchange energy is large, theoretical treatment is not so straightfor- 
ward. Asano et. al have used a recursive Green's function in the lattice model (see Fig. 
20) and calculated pair amplitudes, LDOS and Josephson current. In the actual calcu- 
lation, by changing the exchange potential Vex, they have studied Josephson current in 
S/DN/S, S/DF/S and S/HM/S junctions. At the interface, spin flip scattering is in- 
troduced. On-site site scattering potentials are given randomly. The spin-flip scattering 




Fig. 20. (Color online) (a) A schematic figure of a SFS junction on the tight-binding lattice, (b) 
The density of states for each spin direction. The Josephson junction is of the S/DN/S, S/DF/S, and 
S/ HM/S type for Vex/t ~ 0,1 and 2.5, respectively. Vex is the exchange potential in DF and HM. 
[Reproduced from Fig. 1 of Phys. Rev. Lett, 98,107002,(2007) by Asano et. al] 



at the junction interfaces opens the Josephson coupling via odd-frequency spin-triplet 
Cooper pairs. In the middle of HM, the pairing symmetry is purely OTE state with 
equal-spin-triplet pairing. In Fig. 21(a), the pair amplitudes in the middle of the DF 
and DN are plotted for S/HM/S and S/DN/S junctions, respectively . For S/HM/S 
junction, the pair amplitude is an odd-function of Un and has an OTE symmetry. On 
the other hand, it is an even function for S/DN/S reflecting on the ESE symmetry. 
In Fig. 21(b), corresponding p{u)) in the middle of the DF and DN are plotted for 
S/HM/S and S/DN/S junctions, respectively . The p{uj) in the S/HM/S junction has a 
large peak at the Fermi energy by contrast to that in S/DN/S junctions. Therefore, the 
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odd-frequency pairs can be detected experimentally by using the scanning tunneling 
spectroscopy. ^^^^ 




Fig. 21. (Color online) (a) w„ dependences of pair amplitudes in the middle of HM(half metal) and 
in a DN (diffusive normal metal), (b) Corresponding local density of state p{e) [Reproduced from Fig. 
4 of Phys. Rev. Lett, 98,107002,(2007) by Asano et. al.] 



3. Topology and bulk-edge correspondence 

In §2, we discussed various properties of superconductors in terms of symmetry. In 
this section, we argue topological properties of superconductors. 

3. 1 Quantum Hall effect and TKNN number 

Before going to discuss topological properties of superconductors, we would like to 
see more details of the integer quantum Hall states, which are a representative example 
of topological order. 

The quantum Hall states are realized in two-dimensional electron systems under 
uniform magnetic field perpendicular to the plane. In a commensurate periodic potential 
of the crystal field, each level n is described by the Bloch wave function ipn,k with the 
crystal momentum k in the magnetic Brillouin zone (BZ). The integer quantum Hall 
effects occur when the Fermi energy is located in a band gap, and all the levels below 
it are fully occupied. 

To see the topological characterization, consider the first BZ with — tt < k^, ky < it. 
Because of the periodicity, we can identify the edges of the 1st BZ, i.e., kx = — tt 
{ky = — 7r) and k^ = tt {ky = it) are equivalent. Therefore, the 1st BZ can be regarded 
as the torus T^, and the occupied Bloch wave functions tpn,k defines the mapping from 
this to the U{1) phase of the wave functions. The topological index characterizing 
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IQH state 



superconducting state 



bulk 



gapped (Landau level) 



gapped (Cooper pair) 



edge 



chiral gapless edge state 



gapless ABS 



Table VI. Similarity between integer quantum Hall state and superconducting state. 



this mapping is so called the Chern number (or the TKNN integer) defined by^^"^'"^^^^ 



Here, Ai{k) is the "vector potential" of the C/(l) phase in the momentum space defined 
as 



in terms of the periodic part \un.k) of the Bloch wave function, and Bz — dk^Ay — dkyAx 
is the "magnetic fiux density". It can been seen that Ci is an integer corresponding to 
the winding number for the mapping from to the phase of the Bloch wave function. 
Prom the Kubo formula, the Hall conductance an is found to satisfy eq.(l). Thus 
the quantization of the Hall conductance for IQHE can be naturally explained as the 
quantization of the Chern number itself. 

Recently, the topological characterization has been generalized to superconductors 
in the name of "topological superconductor". The key observation is that there 
exists a close similarity between the quantum Hall states and superconducting states, 
summarized in Table VI: Both states are gapped in the bulk, so we need a finite energy 
to create bulk excitations. On the other hand, we may have gapless states on their 
boundaries, i.e. chiral edge state in quantum Hall states and ABS in superconducting 
state. The analogy between them enables us to apply the topological idea of the quantum 
Hall states to superconductors, while they show originally different physical phenomena 
from each other. 

As discussed in §1.2, for quantum Hall states, there is a relation between the chiral 
gapless edge states and the bulk Chern number Ci. So one can naturally expect that 
a similar bulk-edge correspondence also holds for superconductors. In the following 
sections, we will show that this is the case, and gapless ABSs can be related to bulk 
topological numbers of superconductors. 




(48) 




(49) 



En<EF 
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3.2 Bulk- edge correspondence 

Here we will generalize the idea of the bulk-edge correspondence, or bulk-boundary 
correspondence.^®"' ^^^^ Except for quantum Hall states, we may not use the Hall conduc- 
tance to discuss the bulk-edge correspondence, however, we will show in the following 
that the bulk-edge correspondence can be obtained from a general argument of topology. 

First, we discuss how a bulk-topological number can be defined generally. We sup- 
pose a band description of theory, and the bulk wave functions of electrons or quasi- 
particles are given by Bloch wave functions \un{k)) in the first BZ. As is illustrated in 
Fig. 22, we assume that the Fermi energy is located in a band gap so the system is a 
band insulator in the bulk. 



e(k) 




k 



Fig. 22. Schematic picture for the spectrum of a band insulator. The Fermi energy is located in the 
band gap. All the states below the Fermi energy is fully occupied. [Reproduced from Fig. 2 of Bussei 
Kenkyu 94 311 (2010) by Sato.] 

In this situation, we can introduce a bulk topological number by using occupied 
Bloch wave functions as follows: Let us first note that a Bloch wave function generally 
maps a point k in the momentum space to a point in the Hilbert space. Thus, using 
the occupied Bloch wave functions, in which all the momentum space is filled with 
electrons or quasiparticles, we can map the whole of the first BZ to the Hilbert space. 
Then, depending on a type of wave functions, the image of the BZ may "wind the 
Hilbert space" in some way so that it can not be deformed into a point smoothly. If 
this happens, a bulk topological number is defined as a "winding number" of the image 
of the BZ. 

The above definition of topological number is rather abstract, however it is sufficient 
to derive a general property of bulk topological numbers. The definition infers that the 
topological number takes only discrete integer values since it counts a "winding number" 
of the image of the BZ. However, at the same time, the topological number could change 
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only continuously since it is defined by using the wave functions which can change only 
continuously. From the consistency, we can conclude that the topological number cannot 
change actually as long as no singularity arises. 

Then a question is when a singularity arises. The answer is when the bulk gap closes. 
If the bulk gap closes, we have a gap closing point P. (See Fig. 23.) Since we cannot 
distinguish an occupied state from an empty state at P, the concept of occupied states 
becomes ill-defined, so is the topological number itself. Thus the topological number 
may change discontinuously in this case. 



£(k) 



Ef 




k 



Fig. 23. When the bulk gap closes, a gap closing point P appears. We cannot distinguish an occupied 
sate from an empty state at P. [Reproduced from Fig. 4 of Bussei Kenkyu 94 311 (2010) by Sato.] 

From the above properties of topological numbers, we can derive a generalized bulk- 
edge correspondence. Consider an interface between two different band insulators illus- 
trated in Fig. 24(a). We suppose that a bulk topological number u is defined in both 
sides of insulators, and u ^ in the left hand side while u = in the right hand side. 
Let us see what happens to the topological number when we go from the left to the 
right. We immediately find that the nonzero value of the topological number changes 
abruptly to zero near the interface at Xe- Thus, the above discussion of bulk topological 
numbers implies that the gap of the system must close near the interface. In other words, 
we have a gapless state near the interface. Regarding the vacuum as the topologically 
trivial insulator in right hand side, we have a generalized bulk-edge correspondence: On 
a boundary ( edge or surface ) of an insulator with a nonzero bulk topological number, 
there exist gapless states corresponding to the topological number. 

While the above argument is rather intuitive, we can make it rigorous by intro- 
ducing the confining potential separating the topological phases. ^^^^ The quantitative 
arguments confirm the robustness of the bulk-edge correspondence. 

In the following sections, we will see how the bulk-edge correspondence works for 
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Fig. 24. (Color online) (a) An interface between the insulator with a nonzero topological number 
1/ 7^ and that with = Q. (b) Gap AE of the system. When the topological number v changes at 
X = Xc, the gap closes. [Reproduced from Fig. 5 of Bussei Kenkyu 94 311 (2010) by Sato.] 



superconductors. In the ground state of a fully gapped superconductor, the negative 
energy states are fully occupied while the positive energy states are empty. So we can 
regards a superconductor as a kind of "insulator". This identification enables us to 
apply the bulk-edge correspondence to superconductors. 

3. 3 Topology of Andreev bound state with flat dispersion 

First, we consider a class of superconductors in which the gap function is a single 
component real function. It includes conventional s-wave superconductors, px- 

wave (or p^-wave) superconductors, high Tc cuprates, and so on. When the spin-orbit 
interaction is negligible, and the Cooper pairs preserve the time-reversal invariance and 
a spin in a certain direction, say Sz, this class of superconductors is realized. 

The Hamiltonian is given by 

^ = E(4t>c-H)^(fc) ( I (50) 

k \ ^-ki I 



where 'H(fc) is a 2 x 2 Bogoliubov de Gennes (BdG) Hamiltonian 

e{k) A{k) 
A{k) -6{k) 

with 



_ J = fPi-k) for spin-singlet ^^^^ 

dz{k) = —dz{—k) for spin-triplet 

Here Cfco- (c^o-) annihilation (creation) operator of electron with momentum k 

and spin a, and e{k) the energy dispersion in the normal state. {s{k) = e{—k).) The 
time-reversal invariance implies that the gap function A{k) can be chosen to be real. 
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Diagonalizing the 2x2 BdG Hamiltonian 'H(fe), we find that the quasiparticle spectrum 



E{k) is given by E{k) — ±^/e{ky~+~A(k)^. The gap of the system closes when the 
following condition is satisfied, 

e{k) = 0, A(fe) = 0. (53) 

The negative energy state is fully occupied in the ground state. 

We will identify the Hilbert space of this model. Since the BdG Hamiltonian is a 

2x2 real symmetric matrix, the occupied state \u{k)) is a two component real vector 
with unit norm, which is given by 

, , / cosQ;(fe) \ ^ ^ 

Hk)) =[ , (54) 

\ sm a[k) I 

with an angle variable a{k). Because of an sign ambiguity in the occupied state, 
\u{k)) —\u{k)), the angle variable a{k) should be identified with a{k) + n, 

a{k) ~ a{k) + n. (55) 

In other words, the state with a{k) is physically the same as that with a{k) + vr. Thus, 
instead of the eigenstate (54) itself, the physical state should be rather characterized 
by the unit vector, 

cos 2a{k) 
sin 2a{k) 

The Hilbert space of the model is a one-dimensional sphere parameterized by the 

vector (56). 

Now let us see how to define a topological number in this class of superconduc- 
tors. We first note a subtlety of the system. As was discussed in §3.2, a fully gapped 
system is needed to obtain a well-defined bulk topological number. However, for super- 
conducting states described by eq.(51), we often have a nodal superconductor such as 
two-dimensional p^-wave superconductor or a high Tc cuprate. This is because the gap 
closing condition (53) is rather easily met in the two- or three-dimensional momentum 
space. To resolve this problem, we regard the momenta in certain directions as param- 
eters of the system. Fixing them to certain values, we effectively have a "fully gapped 
one-dimensional system" , in which a bulk topological number can be defined. For con- 
crctcncss, in the following arguments, we fix the momenta in the y and ^-directions, 
and consider the system as a one-dimensional system extended in the x-direction. (See 



(56) 
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Fig. 25.) 



-7^ 




Fig. 25. (Color online) Two-dimensional p^^-wave superconductor. Since nodes exist at A^i and A''2, 
the system is not gapful in the whole of the BZ. However, if we restrict ourself to the momentum space 
with a fixed ky = ky, we obtain a fully gapped one-dimensional system. [Reproduced from Fig. 6 of 
Bussei Kenkyu 94 311 (2010) by Sato.] 



In the one-dimensional gapped system, the BZ is — tt < fc^. < tt. It is essentially a 
one-dimensional sphere (circle) since fca; = vr is identified with = — vr. Thus the 
occupied state \u{k)) maps of the Brillouin zone to of the Hilbert space. Counting 
how many times the image of the BZ winds around the Hilbert space, we can define 
a topological number. Every time the image winds around the Hilbert space, 2a{k) in 
eq.(56) changes by 27r. So the winding number wid is given by 

1 r 

widiky,k^) = —J dk^dk,{2a{k)). (57) 
Using 9{k) defined by 

e{k) 



cos 6{k) 



v/£(fc)2-f- A(fc)2' 



sin^(fc) = , ^^^^ =, (58) 

y/e{ky + A{ky' 

we can show that the occupied state \u{k)) is given by 

/«-.)/2|\ 
\ sinlW*) - ,r)/21 ) 

Thus, we obtain a{k) = 9{k)/2 — 7r/2. From this relation, the topological number wid 
can be rewritten in terms of the parameters of BdG Hamiltonian, 

1 

wid = ^ I dk^e''^ma{k)dk,mb{k), (60) 
where ma{k) is given by 

e{k) 



mi{k) 



^e{ky + A{ky' 
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Then, the above integral is evaluated as a simple sum^^°'^^^'^^^) 

«^id = -^ sgn[A(fc)]-sgn[afc,£(fc)], (62) 

kx;e{k)=0 

where the summation is taken for kx with e{k) — 0. Prom the bulk-edge correspondence, 
there exist the gapless states on the boundary when wid is nonzero. 

We note here that the resultant ABS has flat dispersion. This is because the topo- 
logical number Wid{ky, kz) is nonzero in a finite region of {ky, kz) since it cannot change 
unless the integration path intersects a gap node. From the bulk-edge correspondence, 
this implies that the zero energy state also exists in a finite region of {ky, kz)- In other 
words, the ABS corresponding to nonzero wid has a flat dispersion. 

For the flat dispersion ABS, the bulk-edge correspondence is nicely summarized in 
the form of the index theorem. ^^^'^^^^ Since the BdG Hamiltonian (51) has the so called 
chiral symmetry 

{n{k),ay} = Q, (63) 

it can be shown that the ZEABS is an eigenstate of the chirahty operator ay. Then, 

denoting the number of the zero energy states with the eigenvalue cr^ = ±1 as we 
can relate the index n'^^ — ng ■* to the winding number wid in the form of the index 
theorem, 

wm = - nj,-^ (64) 

or 

wid ^ nl)~^ - n'^\ (65) 

where eq.(64) (eq.(65)) holds for the ABS on the surface of the semi-infinite supercon- 
ductor on X > {x < 0)}^^^ 

Now we will see that the bulk-edge correspondence reproduces the criterion of the 
ZEABS proposed previously. ^^'^'^^^^ In the case where the topology of the Fermi surface 
is simple as illustrated in Fig. 26, it has been known that if the gap function satisfies 

A{kx, ky, kz)A{-kx, ky, kz) < 0, (66) 

then a ZEABS exists on the boundary perpendicular to x-direction.^^'^^^ In other words, 
a sign change of the gap function with respect to — )■ —k^ implies the existence of the 
ZEABS on a surface perpendicular to the a;-direction. The bulk-edge correspondence 
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reproduces this result exactly: Equation (62) leads to 

wid = ~^[sgn[dk,e{-kl,ky,k,)]sgn[A{-kl,ky,k,)] 

+ sgn[dkAkl K, fc,)]sgn[A(fcO, ky, k,)]] , (67) 

where {±k^,ky,kz) denotes the intersection points between the integral path of 
wid{ky) and the Fermi surface. See Fig. 26. Noticing that sgn[(9fc^£:(fc°, /cy, /c^)] = 
—sgn[dk^e{—k'^, ky, k^)], we can rewrite this as 

wid = -hgn[dk,e{-kl,ky,k,)] 

X [sgn[A(-A;0, ky, k,)] - sgn[A(A;0, ky, k,)]] . (68) 

Thus the topological number Wi^ becomes nonzero only when the gap function satisfies 
eq.(66), which means that the bulk-edge correspondence reproduces the previous one 
in this particular simple case. 




Fig. 26. (Color online) Fermi surfaces with simple topology in (a) quasi-one-dimensional system, 
(b) quasi-two-dimensional one, and (c) three dimensional one. The thick blue lines denote the integral 
path of wid- For simplicity, we illustrate the integral path only in (a) and (b). For each case, the 
integral path gets across the Fermi surface only twice at kx = ±^2- (^) (b), P and Pf denote 
the intersection point {—k^,ky) and (k^jky), respectively. [Reproduced from Fig. 4 of arXiv:1102.1322 
to appear in Phys. Rev B by Sato et al.] 

It should be noted here that the bulk-edge correspondence does not merely reproduce 

the known criterion, but is more informative. It is also applicable to more complicated 

cases in which the previous criterion does not work. Such examples were presented in 
refs.160) g^nd 161) 

As a concrete example, we first consider the two-dimensional dxy-^wave superconduc- 
tor where e{k) and A(fc) in eq.(51) are given by 

e{k) = ^-^^, A(fc) = Ao^. (69) 
Here Aq is a positive constant. From eq.(62), the topological number wid{ky) is evaluated 
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as 

{-1, for < ky < kp 
1, for > ky> -kp , (70) 
0, for \ky\ > k-p 

where kp — ^/^/rnjl is the Fermi momentum. Thus the index theorems (64) and (65) 
imply the existence of the corresponding ZEABS. 

The ZEABS is obtained by solving the BdG equation directly For the semi-infinite 
dxy superconductor on x > with the boundary condition \u{x — 0,ky)) — 0, the 
ZEABS on X = is given by^^) 



uo{x))^C\ e^*^^^sin(A;^x)e-^/^ (71) 




where C is a normalization constant, \ky\ < kp, kx = ykp — ky and = mAQky/kp. 
Since the ABS is an eigenstate of ay with eigenvalue ay = —1 {ay = 1) for < ky < kp 
{0 > ky > —kp), it is found that Uq^^ = and tIq ^ = 1 for < ky < kp {n^Q^^ = 1 and 
no~^ = for > /Cy > -kp). When \ky\ > kp, no ZEABS is found, thus n'^^ = Wq"^ = 0. 
As summarized in Table VII (a), the index theorem (64) holds. 

On the other hand, for the semi-infinite dxy superconductor on x < 0, the ZEABS 
on the surface at x = is given by 

\uo{x))^C\ ^ I e*'=^^sin(it^a;)e^/«. (72) 
\ isgnky I 

Thus, the index theorem (65) holds in this case. (See also Table VII (b).) 

Now consider the two-dimensional p^^-wave superconductor. The gap function is 
A(fe) = Aokx/k with k — Vk^ and e{k) is the same as that in eq.(69). For the pj;-wave 
superconductor, we have 





for 






for 





wid{ky) = < , (73) 

I 0, for \ky\ > kp 

from eq.(62). Correspondingly, if \ky\ < kp, we obtain the following ZEABS on x = 



\u 



x,ky))^C{ \ \ e^'^^y sm{kxx)e-''/^^ , (74) 



—I 
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(a) dajy-wavc superconductor on x > 



ky 




4-) 


„(+)_„(-) 


Wld{ky) 


<ky <kF 
> ky > —kp 
\ky\ > kp 



1 



1 





-1 
1 



-1 
1 




(b) d^y-wave superconductor on x < 


ky 








Wldiky) 


0<ky <kF 

0> ky> -fcp 
\ky \ > fcp 


1 






1 




1 

-1 




-1 
1 





Table VII. The number riQ^-* of the zero energy ABSs with the ay = ±1 for (a) the semi-infinite 
dxy-'wavc superconductor on a; > and (b) that on x < 0. For comparison, we also show the topological 
number wid{ky) given in eq.(70). The index theorem (64) and (65) hold in (a) and (b), respectively. 



for the semi-infinite p-^-wave superconductor on x > 0, and 

\u{x,ky)) = C 1^ ^ j e*^^^sin(A;^x)e^/^^ (75) 

for the semi-infinite Px-wave superconductor on x < 0. Here C is a normahzation 
constant, kx — ■^k^ ~ ky and Cp^ — inAo/kp. It is also found that these solutions are 
the eigenstates of ay with the eigenvalue ay — —1 and ay — 1, respectively. Thus n'^^ 
and tt-q"^ are summarized as Table VIII (a) and (b). We confirm the relations (64) and 
(65), respectively again. 

Finally, we would like to point out the relevance between the topological structures 
of the ABS and the odd-frequency pairing. The surface odd-frequency pair amplitude 
discussed in §2.2.2 is rewritten as 

- fe^ ^^^^ 

for spin-singlet dxy-^a^Q superconductor and 

h = (77) 



for spin-triplet p-wave superconductor by cos6 and ky = kpsinO. As 

well as the wave function of ABS of spin-singlet dxy and spin-triplet p^-wave super- 
conductor derived in this section, the factor sgnky exists only for spin-singlet dxy-wetve 
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(a) px-wdvc superconductor on a; > 



ky 




"-0 


"■0 "0 


Wld{ky) 


\ky \ < k-p 
\ky\ > kp 






1 



-1 



-1 




(b) Px-'wslvc superconductor on x < 


ky 




"-0 


J+) _ J.-) 


Wldiky) 


\ky \ < fcp 
\ky \ > fcp 


1 







1 



-1 





Table VIII. The number of the zero energy ABSs with the ay = ±1 for (a) the semi-infinite 

Px-wave superconductor on a; > and (b) that on x < 0, respectively. For comparison, we also show 
the topological number Wid{ky) given in eq.(73). The index theorem (64) and (65) hold in (a) and (b), 
respectively. 



case. This factor decides the difference of the parity of induced Cooper pair. Next, we 
consider a two-dimensional semi-infinite superconductor in x < 0. The corresponding 
pair amphtude at surface {x — 0) is given by 

? isgn{ky)AQ \K\\ky\ 

~ ^ ^ 

for spin-singlet dxy-wave superconductor and 

for spin-triplet p^-wave one, respectively. Comparing eq. (76) [(77)] with (78) [(79)], it 
is evident that the difference between them is the presence of — sign. The present — 
sign exactly corresponds to the different values of and rio '* between case (a) and 
case (b) in Tables VII and VIII. 

3.4 Time-reversal breaking superconductors and Majorana fermion 

Next, consider time-reversal symmetry breaking superconducting states. The anal- 
ogy between the quantum Hall states and superconducting states is direct in this 
case.^^^'^^^'^^^'^^^^ The simplest one is a superconductor with a single complex gap func- 
tion. As well as the previous subsection, it is realized when the spin-orbit interaction is 
negligible, and the Cooper pairs preserve Sz, but break the time-reversal invariance in 
this time. Also, it includes a spinless superconductor where the Cooper pairs are formed 
by fully spin polarized electrons. We consider the latter case first. 
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For the spinless superconductor, the Hamiltonian is given by 

Cfct 



(80) 



with 



KWH '^'^ |. (81) 

' A*{k) -€{k) ' 

Here we assume that the electron is fully polarized with up spin. From the Fermi 
statistics, A(fc) is an odd function of k. Thus the superconducting state is spin-triplet. 

The above BdG Hamiltonian H{k) is a 2 x 2 hermitian matrix, so the occupied state 
\u{k)) is a two-dimensional complex vector with unit norm, 

( cosa(fe)e-''^('=) \ 
^''^'''^^ y sina(fc)e-^W ) ' ^^^^ 

where a{k), f3{k) and 7(fe) are angle variables with < a{k) < 7r/2, < /3(fe) < 27r 
and < j{k) < 27r. Using a phase ambiguity of the eigenstate, 

\u{k)) ^e"^^^^\uik)), (83) 

we can set 7(fc) = in eq.(82). To identify the Hilbert space, we calculate the expecta- 
tion value of the Pauli matrices Uj, 

i{u{k)\a^\u{k)), {u{k)\ay\u{k)), {u{k)\a,\u{k))) 

= (sin2a(fc)cos/3(fc),sin2a(fc)sin/3(fe),cos2a(A;)). (84) 

The vector (84) parameterizes the two-dimensional sphere S'^. Thus we have a one- 
to one correspondence between the occupied state \u{k)) and a point on S^. In other 
words, the Hilbert space in this model is S'^. 

In two dimensions, the BZ is topologically equivalent to the two-dimensional torus 
T^. Thus the occupied state \u{k)) maps of the BZ into 5*^ of the Hilbert space. 
Since the surface element of of the Hilbert space is 

sin2Q;d(2Q;)d^, (85) 

the winding number of the image of the BZ is evaluated as 

W2d = ^j J dhdkySm[2a{k)]€'^dkMk)dk,P{k). 

(86) 
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Introducing 9{k) and (p{k) as 

sik) 



cos9{k) 



V£(fc)2 + |A(fc)P' 



sing(fc)e^^W= , ^^^^ =, (87) 
Ve(/e)2+|A(fe)P' 

the occupied state is given by 

/cos«,W/2-,/2)e<^<'.\ 
\^ sin(e(fc)/2-7r/2) y 

which leads to 

a{k) =e{k)/2-7r/2, (3{k) = -ip{k). (89) 
Prom this, the winding number is recast into 

dk:cdkye'h'''"'ma{k)dkimb{k)dkjmc{k), 



"Stt 



where ma{k) is given by 

ReA(fc) 



mi(fc) 
m2{k) 



ImA(fe) 

V£(fe)2+iA(fe)F' 



(90) 



m^ik) = , ^^^^ ■ (91) 

V£(fc)2 + |A(fc)P 

Finally, this can be evaluated as a simple sum^^°' ^^^^ 

W2d^~ J2 sgn[s{k)]sgn[detdkAj{k)], (92) 

A(fc)=0 

with Ai(fc) = ReA(fc) and A.2{k) = ImA(fc).^^®) Here the summation is taken for 
{k^, ky) with A(fc) = 0. When W2d is not zero, we have a chiral edge state on the 
boundary. 

In Fig 27, we illustrate the quasiparticle spectra of spinless chiral p-wave supercon- 
ductors. Here we have used the lattice model with 

e{k) = —2tx cos kx — 2ty cos ky — II, 

A(fc) = Ao(sin /ca; + isin /cj^), (93) 
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and the quasiparticle spectra are calculated under the open boundary condition at 
X = and x = L in the x-direction. The resultant winding number W2d depends on 
the topology of the Fermi surface, and correspondingly, the gapless edge modes also 
depend on the topology of the Fermi surface. For a two-dimensional Fermi surface, we 
find that W2d = 1, so we have a single chiral edge mode on each side of boundary. 
The relation between the Fermi surface topology and the winding number W2d (and the 
corresponding edge state) is a general property of spin-triplet superconductors, ^^'^^ as 
will be explained in §3.7. 




Fig. 27. (Color online) Edge states in spinless chiral p + ip-wave superconductors given by eq.(93). 
We take (a) = 0.2, ty = 1, ^ = -1, Aq = 0.5, (b) = 1, ty = 1, = -1, Aq = 0.5, and (c) 
tx — ty = 0-2, = —1, Aq = 0.5. (Top panels) Fermi surfaces in the normal states. (Middle panels) 
Quasiparticle spectra of the corresponding superconducting state with edges at x = and x = L. The 
gapless modes are localized on the edges. The winding number W2d is (a) W2d = 0, (b) W2(i — 1, and 
(c) W2d = 0, respectively, (bottom panels) Schematic illustration of the corresponding edge state. The 
arrows indicate the direction of the group velocity of the edge modes. 



As well as the edge state, we can apply the bulk-edge correspondence to the vortex 
in superconductors. A vortex in a superconductor can be regarded as a hole in 

the bulk. Then, applying the bulk-edge correspondence to the edge of the hole, we find 
that there is a zero energy state in the vortex core when W2a 0. Some zero modes 
become non-zero mode when deforming the hole into a real vortex. When W2d is odd, 
however, at least a single zero mode survives since the zero modes become massive in 
a pair due to the particle-hole symmetry. In particular, there is a single zero mode in a 
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vortex when 'W2d — ±1.^^^"^^^^ 

Here we would like to mention that the analogy between the quantum Hall states and 
superconducting states is direct in this case. The quantum Hall states are topologically 
characterized by nonzero Chern number. For comparison, let us evaluate the Chern 
number in this model. By using the occupied state (82), the Chern number is given as 

= ^ /" f dk^dkye'^dk,Aj{k). (94) 

27r J_^J_^ 

where Ai{k) = i{u{k)\dkiU{k)) . From a straightforward calculation, we obtain the rela- 
tion 

e'^dk^Aj = ~e'^e'''''madk,mbdk,mc. (95) 

This implies that winding number W2d is nothing but the Chern number, W2d = Ci. 

In spite of the analogy above, there are differences between them. The U{1) elec- 
tromagnet gauge symmetry is spontaneously broken in the superconductor while it is 
not in the quantum Hall states. Thus no quantized Hall current is carried by the edge 
state in the former. Furthermore, the quasiparticle in the superconductor 

- I I , (96) 



satisfies the additional condition. 



*Lfc = r^'fc, r= 1^ ° J j , (97) 

The condition (97) is called the Majorana condition. It means that the quasiparticle 
is essentially the same as its antiparticle ^L^- As a result, the edge state in the 

superconductor is a 1 -|- 1 dimensional Majorana chiral fermion, not a chiral fermion in 

the quantum Hall state. 

An important consequence of the Majorana condition is that it gives rise to the 

non-Abehan statistics of the vortices. ^^^'^^^^ As shown in the above, for the spinless 

chiral p-wave superconductor, there is a single zero mode 70 in the vortex core. The 

zero mode satisfies the anticommutation relation 

{7^,70} = 1- (98) 
At the same time, the Majorana condition reads 

70 = 7^- (99) 
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Thus 7o is neither the annihilation operator nor the creation operator. To obtain the 
well-defined creation operator, we need a pair of vortices. Using the zero mode 7o*^ of 
the vortex i {i — 1,2), we can construct the creation operator ^'^ as 

y = ^° , (100) 

which satisfies 

7 7^7^ {7,7^} = 1- (101) 

Since the vortex 1 and 2 are separated from each other, the creation operator is defined 
non-locally. This non-locality changes the statistics of vortex drastically. Indeed, the 
vortices obey the non-Abelian statistics. 

To sec the non-Abelian statistics of the vortices, consider the vortex 1 and 2 il- 
lustrated in Fig. 28. Then, the vortex 3 encircles the vortex 2 and it goes far away. In 
classical theory, the final configuration is completely the same as the initial one. In 
quantum theory, however, the final state can be completely different from the initial 
one. 

Suppose that the initial state |0) is annihilated by 7, 7|0) = 0. When the vortex 3 

encircles the vortex 2, the vortex 2 also encircles the vortex 3 in the rest frame of the 

(2) 

vortex 2. Thus the zero mode 7q in the vortex 2 gets the Aharanov-Bohm phase in 
this process. If an electron (a hole) moves around the vortex 3, the Aharanov-Bohm 

phase is e*^*" = — 1 (e"**^*" = —1), where $0 is the unit magnetic fiux of the vortex, 

(2) 

^0 = Ti'/e. Thus the zero mode 7q gets the same factor, 

7?^ ^ -7? \ (102) 
since it is a superposition of an electron and a hole. This means that 7 changes as 

(1) , ■ (2) (1) • (2) 

7^ ^° ^y = ^o -;^o . (103) 

^2 ^2 

Therefore, the final state is annihilated by 7''', 7^|1) — 0, not by 7. The final state |1) is 

completely different from the initial one |0). Indeed, there is no overlap between them, 

i.e. (0|1) = 0. 

Note that such an exotic phenomenon never occurs if the vortices is boson, fermion, 
or Abelian anyon. For these cases, the encircling process considered here results in a 
overall phase at most, so the final state is essentially the same as the initial one. In 

other words, the above phenomenon is a result of the non-Abelian statistics of vor- 
tices. A particle obeying the non-Abelian statistics is called non-Abelian anyon, and a 
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a) b) c) 




Fig. 28. Movement of the vortex 3 around the vortex 2. 

topological phase with non-Abelian anyons is called non-Abelian topological phase. 

Up to now, we have focused on the spinless superconductor. Now consider the spinful 
time-reversal breaking superconductor. An example is a chiral p-wave superconductor 
in Sr2Ru04 or in a ^He A-phase superfluid thin film. When the spin-orbit interaction 
is negligible and the Cooper pairs preserve a spin in a certain direction, say Sz, the 
Hamiltonian is given by 

'^ = E(4t'C-H)^(fc) f I' w 

k \ c^ki I 

with 'H(fc) in eq.(81). The gap function A(fc) is given by eq.(52). Since ^{k) is the same 
as that for the spinless superconductor, the Hilbert space is also the same, that is S'^. 
Thus the topological number characterizing the state is the same winding number W2d- 
The difference between the spinless superconductor and the spinful one is the Majorana 
condition. For the spinful superconductor, the quasiparticle operator \E'fe 




(105) 



in eq.(104) does not satisfy the Majorana condition (97). Nevertheless, this does not 
imply that the Majorana fermion cannot appear in the spinful superconductor. Indeed, 
in a special kind of vortex called half-quantum vortex, the Majorana fermion appears. 

The half-quantum vortex is realized in spinful spin-triplet superconductors, where 
the spin components of the gap function can rotate around vortex. Using this degrees 
of freedom, the half-quantum vortex has the following configuration, 

d{k) = A{k, r)e^^/2 (cos 6/2, sin 6/2, 0) , (106) 

where d{k) is the d-vector of the spin-triplet gap function, r the distance from the vortex 
core, and 6 the angle around the vortex. Since this configuration does not preserve 5*^, 
its BdG Hamiltonian is not given in the form of eq.(104). Instead, we can treat the 
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system as two copies of spinless superconductors. Indeed, the matrix representation of 
the gap function leads to 



so there is no mixing between the different spin sectors. Then, each spin sector is the 
spinless superconductor described by eq.(80). (For the down spin sector, and 
should be replaced by Ck^ and c^, respectively.) Equation (107) implies that only the 
down-spin sector has a vortex. Therefore, we can regard as the half quantum vortex as 

the vortex in spinless superconductor with down spin. This means that if W2d is odd, 
there is a Majorana zero mode in the half-quantum vortex. Thus the half-quantum 
vortex can obey the non-Abclian statistics. Recently, an experimental evidence of the 
half-quantum vortex in Sr2Ru04 has been reported. ^^^^ Also, a half-quantum vortex is 
expected to exist in a thin film of ^He A phase. ^'^^'^'^^^ 

3.5 Helical edge modes in TI/QSHE 

In this subsection, we consider the band insulators without the superconductivity 
from the viewpoint of the topology. For time-reversal invariant systems, new topological 
phases arise when the spin-orbit interaction is strong enough. Here we review briefly such 
states, the quantum spin Hall states^^^"^^^^ and its three dimensional generalization, 
topological insulators. ^^^"^^^^ 

The quantum spin Hall state is easily understood as a pair of the integer quantum 
Hall states if is a good quantum number. Consider the following two-dimensional 
system preserving S^, 



where T-L+{k) {l-L-{k)) is the Hamiltonian for up-spin (down-spin) electron. Here we 
suppose that the spin-orbit interaction works as a kind of magnetic field, and each 
spin sector realizes an integer quantum Hall state with the Hall conductance cr^^ = 
—C^^e^/h (cr = ±). Due to the time-reversal invariance, the Chern numbers of the 
two spin sectors are opposite from each other, = —C^\ Thus, the spin Hall 

conductance (JsH = o"]^^ — o'h quantized, while the total Hall conductance, an = 
a^jj^^ + a\j \ is zero. Correspondingly, there exist C['^^ {C[ ^) gapless edge states with 
up-spin (down-spin) carrying the spin Hall conductance. 




(107) 




(108) 
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Now consider a time-reversal invariant perturbation which does not preserve Sz- 
The spin indices a — ± turns into that of pseudo-spin a, and the Chern numbers 
Ci^^ and C[~^ are replaced by C["^ and C[^\ Then, due to the spin mixing, most of the 
gapless edge states become massive. Nevertheless, if C[°'\= —C[^^) is an odd number, 
at least one pair of gapless edge states survives: Because of the time-reversal invariance, 
the gapless states form a set of Kramers doublets. Since there is no mixing in the same 
Kramers doublet, the Kramers doublets become massive in a pair by the mixing with a 
different Kramers doublet. Thus, if there are an odd number of the Kramers doublets, 
which is reahzed when C[°'^ is odd, at least one Kramers doublet of the gapless edge 
states remains massless without a partner. 

As indicated by the above argument, for a generic time-reversal invariant system, 
where is not a good quantum number, the quantum spin Hall state is characterized 
by the parity of C|"\ (^_i)'^2dTi = (^_i^c[ ^ ^ which is called the Z2 invariant in the 
literature. Then the resultant gapless edge state is called the helical edge mode. 

The topological insulator is a three dimensional generalization of the quantum spin 
Hall state. ^'^^"^'^''^ It is characterized by the four Z2 invariants {uq; vi, 1/2, vs), which can 
take the value or 1. When, i/q = 1, the system is called strong topological insulator 
which supports an odd number of 2D Dirac fermions on the surface. When some of ui, 
1/2 and 1/3 are 1, while uq — 0, the system is classified into weak topological insulator. In 
this case, an even number of Dirac fermions are on the surface, which can in principle 
be paired and gapped due to the disorder scattering etc., and are not stable. However, 
the weak topological insulator is still distinct from the conventional insulator with all 
z/'s being zero, as has been evidenced by one-dimensional channels appearing along the 
dislocation. ^^^^ 

3.6 Helical superconductors and helical Majorana modes 

There exist a topological state of superconductors which is analogous to the quan- 
tum spin Hall state. Such a superconductor is called the helical superconductor. ^^^'^^^^ 
A representative example of the helical superconductor is two-dimensional Rashba non- 
centrosymmetric superconductors. ^^^'^^^^ Let us see its basic property. 

The model Hamiltonian of the Rashba superconductor in two dimensions is 



— "Hkin + ^so + ^ 



•pairing ) 
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k,cr,a' 



H 



■pairing 




+ h.c, 



(109) 



where c^^ {cka) is a creation (an annihilation) operator for an electron with momentum 
k = (kxyky), spin a. The energy band dispersion is £k = — 2t(cosA;^ + cos ky) — /j, 
with the hopping parameter t and the chemical potential fi, and the Rashba spin-orbit 
coupling is jC^ = {sinky, — sink^). Because of parity mixing of Cooper pairs, the gap 
function A{k) has both a spin-triplet component d{k) and a spin-singlet one ih{k) at 
the same time, A(fc) — iilj{k)ay + id{k)aay. Due to the strong spin-orbit coupling, 
the spin-triplet component d{k) is ahgned with the Rashba coupling, d{k) — At/^fc.^^^^ 
For the spin-singlet component ip{k), we assume an s-wave pairing, ip{k) — Ag. The 
amplitudes At,s are chosen as real. The Zeeman coupling /ibHz X]fe(clt^''t~^l4,^fe4-) with 
Hz a magnetic field in the z direction has been also introduced for later use. 

Before going to study topological properties of the system, we first examine the 
bulk spectrum of the system. Topological nature of the system changes only when the 
gap of the bulk spectrum closes. The bulk spectrum E{k) of the system is obtained by 
diagonalizing the following 4x4 matrix. 



n{k) 



£k - ha^ + aCk ■ cr iA^o-y + iAtCk ■ cray 
-iAsGy - iAtCkCTy ■ a -Sk + ha^ + aCk ■ cr* 




(110) 



with h = ij,bHz, and we have 



±el + {a' + Al)Cl + h' + Al 




(111) 



The gap of the system closes only when 



4 + {a' + Al)Cl + h' + Al 




(112) 



J. Phys. Soc. Jpn. 



FULL PAPERS 



which is equivalent to 

£feAt A = A,aCk. (113) 

When At 7^ 0, eq.(113) is met either when 

= (^l + ^yA^Cl-AD+h'^O. (114) 

or 

el + Al^h\ A = 0. (115) 

In the absence of the magnetic field, only eqs.(114) can be met and they are rewritten 

in simpler forms, 

el^a^Cl AlCl^Al (116) 

Topological nature of the system does not change unless eq.(114) or eq.(115) (or eq.(116) 
when Hz = 0) is satisfied. 

When Hz = 0, the system is time-reversal invariant, and the topological property is 
characterized by the Z2 invariant like the quantum spin Hall state. Below, we will show 
that if the spin-triplet pairing is stronger than the spin-singlet one, the Z2 number is 
non-trivial. 

To see this, we adiabatically deform the Hamiltonian of the system without gap 
closing. This process does not change the Z2 topological number, since it changes only 
when the gap closes. ^^^'^^^^ Prom eq.(116), it is found that if the spin-triplet amplitude 
A^Ck is larger than the spin-singlet one Ag on the Fermi surface given by Ek — ajCo{k), 
we can take Ag — >■ 0, then a — >■ without gap closing. (If = at one of the 
time-reversal momenta k = (0, 0), (tt, 0), (0, tt), (tt, vr), the gap closes when As = 0. 
However, this undesirable gap closing can be avoided by changing /i or t slightly.) Thus, 
its Z2 number is the same as that of the pure spin-triplet SC with d{k) = A^Ck. 
The resultant system preserves Sz, and the BdG Hamiltonian is decomposed into two 
spinless chiral superconductors in up-spin sector and the down-spin sector, 7^+(fe) and 
H-{k), respectively, 

n±{k) = 

£fc —At (sin kx^i sin ky) 

-At (sin kx T i sin ky) -£k 
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Fig. 29. (Color online) The energy spectra of the 2D noncentrosymmetric superconductor with 
edges in the absence of magnetic field. We take t — 1, fi = —3, a = 0.6. (a) noncentrosymmetric 
superconductor with dominating p-wave paring. At = 0.6 and As — 0.1. (b) noncentrosymmetric 
superconductor with purely s-wave paring. At = and Ag — 0.6. 



(117) 

If the Fermi surface is two-dimensional and electron-like, their Chern numbers is Cj^'' = 
=pl for At > 0. (see Fig.27). Thus the Z2 number is non-trivial. 

From the bulk-edge correspondence, there should exist gapless edges if the spin- 
triplet pairs dominate the superconductivity. In Fig. 29 (a), we show the energy spectrum 
of the 2D noncentrosymmetric superconductor with edges. It is found that there exist 
gapless edge states in the bulk gap. The gapless edges states form a Kramers pair. 

For comparison, we also illustrate the energy spectrum for the 2D noncentrosym- 
metric superconductor with purely s-wave paring in Fig.29 (b). As is seen clearly, no 
edge state is obtained. This is also consistent with the trivial Z2 number of the purely 
s-wave paring. 

We also notice that the helical Majorana gapless edge states are very sensitive to 
the direction of the applied magnetic field. As seen in Fig. 30, they become unstable 
under a small magnetic field in the y-direction, while the gapless edge states are stable 
under a magnetic field in the x- and z-direction. As a result, the magnetic field Hy along 
the edge causes a tiny gap of the order 0{fiBHy) for the edge states. The experimental 
detection of the edge states have been proposed theoretically.^^^'^^^'^®^"^^^-' 

3.6.1 Non-Abelian topological order and Majorana fermion induced by spin-orbit in- 
teraction and Zeeman field 
The spin-orbit interaction enables us to realize Majorana fermion and non-Abelian 
topological order in the presence of the Zeeman field. ^^^'^^^'^^^^ In particular, the Ma- 
jorana fermion can be realized even if the symmetry of the gap function is spin-singlet 
dominant. ^^^"^^^^ Here we see this interesting phase of noncentrosymmetric supercon- 
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Fig. 30. (Color online) The dependence of gapless edge states on the direction of the magnetic field. 
We take t = 1, = -3, a = 0.6, At = 0.6 and = 0.1. a) fiBH^^ = 0.15, ^iBHy = and ^ibH^ = 0. 
b) IjLbHx = 0, jjL-QHy = 0.15 and /ie-ffz = 0. 



ductors. 

An intuitive understanding on this mechanism of non-Abehan topological phase is 
obtained if we perform the dual transformation of the BdG Hamiltonian.^^^'^^^^ For the 
two-dimensional Rashba superconductor (109), the dual Hamiltonian is obtained by the 
following unitary transformation, 

H'^ik) = DH{k)D\ D = ^[ ^ '"^M , (118) 



V2 



which leads to 




:ii9) 



- At£fc ■ cr - haz -iSkCTy - iaCk ■ cray 
iSkC^y + iaCkCTyCr -As + AtCk ■ cr* + ha^^ 

From eq. (119), it is found that the Rashba spin-orbit interaction in the original Hamil- 
tonian is formally transformed into a "p-wave pairing interaction" with the d vector, 
d{k) = —aCki in the dual Hamiltonian l-(P{k). Since 'HP{k) has a nonstandard con- 
stant kinetic term, this does not necessarily mean that the topological properties of 
T-L{k) are the same as those of a usual p-wave superconductor. However, we find that 
that the topological order emerges when h satisfies 

h^>/\l + el^^, (120) 

due to the Rashba spin-orbit interaction. Indeed, one can prove that the Chern number 
Ci = —1 if the condition (120) is niet.^^^'^^^'^^^-* Also, the same condition (120) is 
obtained from the analysis of Majorana zero mode in a vortex. ^^^^ In Fig.31, we illustrate 
the energy bands of our Hamiltonian (109) with edges. This figure clearly shows the 
existence of a chiral gapless state if the Zeeman field h satisfies eq.(120), even when the 
spin-triplet superconductor is absent. 
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Prom the argument in §3.4, it is evident that a vortex supports a Majorana zero 
mode in this case. Indeed, we can demonstrate the existence of the Majorana zero mode 



in §3.4, if there exists a single Majorana fermion zero mode for each vortex, vortices 
obey the non-Abehan statistics. Following ref.,^^^-* we use here the dual Hamiltonian 
H^{k) to solve the BdG equation, then construct a solution in the original Hamiltonian 
1-L{k) by using the duality transformation (119). For simplicity, we assume ek=o — and 
At = 0. Then, if h satisfies |As| < \h\ < 2|As|, the low energy properties are governed 
by quasiparticles on the smaller Fermi surface, which is split from the larger one by 
the spin-orbit interaction. ^^^^ The larger Fermi surface can be neglected to construct 
the zero mode. Thus, we concentrate on fermions with |fe| ~ for which yP{k) is 
decomposed 



unique zero energy solution with a quasiparticle field 7^ = / dr[uoiJj\_ + ?^^+], where 
uo = i(re^'')-i/2g-(/i-A,)r/a^ ^ _.^^^-i0yi/2^-ih-As)r/a^ r^^^ sohition is nomializable 

when eq.(120) is satisfied. This is the Majorana zero energy mode; i.e., ^'^ = 7. Using the 
duality transformation (119), we also found that a vortex in the original Hamiltonian 
has a single Majorana zero mode, which implies that the vortex is a non-Abelian anyon. 

Here note that the condition (120) implies the Zeeman energy larger than the s- 
wave BCS gap Ag. Thus the orbital depairing effect could destroy the superconductivity 
while the Pauli depairing effect is suppressed in the presence of the Rashba spin-orbit 
interaction. If the spin-triplet amplitude is dominant and the s-wave pair ampli- 
tude is small enough, the condition (120) can be met for a weak magnetic field without 
destroying the superconducting gap.^^^^ On the other hand, if the s-wave pair amplitude 
is dominant, one always needs a strong magnetic field larger than the superconducting 
gap. Several schemes to realize the non-Abelian topological phase in this case has been 
proposed, (i) s-wave superfiuids in neutral ultracold fermionic atoms with laser gener- 
ated spin-orbit interaction, ^^^'^^^^ (ii) heterostructure semiconductor device, and 
(iii) heavy fermion systems. ^^^^ 

If we consider dj.2_j^2-wave or dxy-wave pairing instead of the s-wave pairing, the 



explicitly by solving the BdG equation for a single vortex. ^^^'^^^ 197,199) explained 




(121) 



The BdG equations for with a single vortex can be solved, ^^^^ then, we find a 
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Fig. 31. (Color online) Majorana edge modes in an s-wave superconducting state with Rashba spin- 
orbit interaction under the Zeeman magnetic field h in the z-direction. Energy spectra of the system 
with open boundary in the x-direction and the periodic boundary condition in the y-direction. (a) 
h — (b) h = he = \J el^o ^ (c) h > he. When the Zeeman field h is larger than the critical value 
he, the Majorana zero mode appears. [Reproduced from Fig.l of Phys. Rev. Lett. 103, 020401 (2009) 
by Sato et al] 

condition (120) changes as^^^^ 

> 4=0, (122) 

since these pair amphtudes vanish at fc = 0. Thus now the condition (122) is inde- 
pendent of the superconducting gap. This means that the Majorana fermion and the 
non-Abehan topological phase could be realized in a weak magnetic field without de- 
stroying superconductivity for systems with small Sk = 0. In Fig. 32, we illustrate the 
Majorana edge modes in this case. The Majorana fermion state is topologically pro- 
tected in spite of the presence of bulk gapless nodal excitations. Due to the existence of 
bulk nodes, one cannot obtain a well-defined Chern number, but the particle-hole sym- 
metry makes the parity of the Chern number well-defined. ^^^-^ The non-Abelian nodal 
superconductor is realizable in an interface between a centrosymmetric nodal supercon- 
ductor such as high-Tc cuprates and a semiconductor. In such a system, because of the 
considerably large superconducting gap, the experimental detection of Majorana modes 
may be easier. 

3.6.2 Dispersionless Majorana fermion 

For a certain type of noncentrosymmetric superconductors, novel Majorana fermion 
is possible. It preserves the time-reversal symmetry, however, unlike the helical 
Majorana fermions, it consists of a single branch of fiat dispersion. 

The time- reversal invariant Majorana fermion is realized for the dxy+p-wave Rashba 
noncentrosymmetric superconductor. Such a pairing symmetry has been suggested for 
heterointerface LaAlOa/SrTiOs.^"^^ Figure 33 illustrates the quasiparticle spectrum for 
the semi-infinite d^y + p-wave Rashba noncentrosymmetric superconductor on x > 0. 
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Fig. 32. (Color online) Majorana fermions in nodal superconductors, (a) Possible realization scheme 
to realize nodal superconductor with Majorana fermion. Due to the proximity effect and the potential 
gradient at the interface, d-wave superconducting state with the Rashba spin-orbit coupling is realized 
in the interface between the high Tc cuprate and the semiconductor, (a-1) Side view (a-2) Top view, 
(b) Energy spectra of the systems with open boundaries at x = and x = L in the x-dircction and 
the periodic boundary condition in the j/-direction for (b-1) c?a.2_j^2-wave pairing, and (b-2) dxy-wave 
pairing. Majorana gapless edge modes at a; = and x = L are depicted, respectively, in green and red 
curves, (c) Schematic illustration of Majorana edge modes counterpropagating on two opposite edges 
for (c-1) (ij,2_j,2-wave pairing and (c-2) d^y-waMe pairing, [(b) and (c): Reproduced from Fig.l of Phys. 
Rev. Lett. 105, 217001 (2010) by Sato and Fujimoto] 

Due to the Rashba spin-orbit interaction, the Fermi surface is spht into two with the 
Fermi momenta ki and (^i < ^2)- See Fig. 33. There exists a single branch of zero 
energy edge state in the momentum region ki < \ky\ < k2 between the Fermi sur- 
faces. The wave function for the zero energy edge state "^miky) can be written as 
^^^rniky) = {ui{ky),U2{ky),Vi{ky),V2{ky)) with Ui{ky) = vl{-ky), U2{ky) = v*{-ky). 
The Bogohubov quasiparticle creation operator for this state is constructed in the usual 
way as 7 '''(/;; J,) = Ui{ky)c^^{ky)+U2{ky)c]^{ky)+Vi{ky)c^{—ky)+V2{ky)ci{—ky), and satisfies 
the Majorana condition 7^(A;y) = •y^—ky). 

Unlike Majorana fermions discussed in previous sections, the present single Ma- 
jorana bound state at ky is realized with time-reversal symmetry. The time-reversal 
invariant Majorana bound state has the following three characteristics, (a) It has a 
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unique flat dispersion: To be consistent with the time-reversal invariance, the single 
branch of zero mode should be symmetric under ky — > —ky. Therefore, by taking into 
account the particle-hole symmetry as well, the fiat dispersion is required. On the other 
hand, the conventional time-reversal breaking Majorana bound state has a linear dis- 
persion, (b) The spin-orbit coupling is necessary to obtain the time-reversal invariant 
Majorana bound state. Without spin-orbit coupling, the time-reversal invariant Ma- 
jorana bound state vanishes, (c) The time-reversal invariant Majorana bound state is 
topologically stable under small deformations of the Hamiltonian preserving ky. The 
topological stabihty is ensured by the topological invariant l^^^^'^^^) 



W{ky) = i-Im 



dkxdk^ In detq{k) 



(123) 



with q{k) — [sk — i'^{k)\ay -\- [aC^ — id{k)] • auy. The flat dispersion of the Majo- 
rana fermion is terminated in a gap node. The last property is closely related to the 
topological stabihty of line nodes of noncentrosymmetric superconductors. 

It is noted that the appearance of surface flat bands in three-dimensional non- 
centrosymmetric superconductors or topological media have been also discussed re- 
cently ^o^-^n) 



3. 7 Topological phase in odd-parity superconductors 

In §3.4, we have seen that the topological properties of the spinless superconductor 
changes as the topology of the Fermi surface changes. This is a general property of 
spin-triplet odd-parity superconductors. ""^^^'^^^^ 

To see this, consider a general Hamiltonian H for odd-parity superconductors. 



H 



^E(4«,c-.«)i/(/c)(' 



where c^^ (cka) denotes the creation (annihilation) operator of electron with momentum 
k. The suffix a labels other degrees of freedom for electron such as spin, orbital degrees 
of freedom, sub- lattice indices, and so on. £(k) is the Hamiltonian of the electron in 
the normal state. The system in the normal state is supposed to be invariant under the 
inversion Cka — ^ J2a' Paa'C-ka' with = 1, P^S{k)P = S{—k). We also assume that 
the parity of the gap function A(fc) is odd, P^A{k)P* = —A{—k). 
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Fig. 33. (Color online) (Top) Angle- resolved LDOS for dxy + p-w&Ye pairing as a function of ky. 
At 7^ and A^. = . (Bottom) The corresponding Fermi surfaces of the system. The Fermi surfaces 
are split into two by the Rashba spin-orbit interaction. Single branch of zero energy ABS exist in the 
region of fci < \ky \ < k2- 



First, consider the case with the time-reversal symmetry breaking. The BdG Hamil- 
tonian (124) has the following particle-hole symmetry, 

CH(k)C^ = -H*(-k), C = I ° M . (125) 

From this, we can say that if \un{k)) is a quasiparticle state with positive energy 
En{k) > satisfying H{k)\un{k)) = i?„(fc)|u„(fc)), then C|m*(— fc)) is a quasiparti- 
cle state with negative energy —En{—k) < 0. We use a positive (negative) n for \un{k)) 
to represent a positive (negative) energy quasiparticle state, and set 

\u,n{k))=C\ul{-k)). (126) 

Let us define the gauge fields 

At\k) = ^Y,Mk)\^k^u^m. A^"^(fc) =^$^K.(fc)l^fcJwn(fc)), (127) 

n>0 n<0 

where the relation 

/^;\k) = A^r\-k) (128) 
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holds from (126). In addition, we have 

A\-\k) + A\'^\k) =itr[U\k)dk,U{k)] = idk An det U{k), (129) 
with the unitary matrix U{k) given by 

Unm{k) = K\k) (130) 

where u!^{k) is the m-th component of \un{k)) — {••■, w^(fc), li^(fc), • • • 

Now consider the time-reversal invariant path Cij in the BZ, which passes through 
the time- reversal invariant momenta Fj and Tj. See Fig. 34. Note that Cij is a closed 
path because of the periodicity of the BZ. The periodicity of the BZ imphes that the 
integration of A^ \k) along is quantized as, 

j dkiA^r\k) = ]-<[> dki[A^r\k)+A^^\k) 

J Cij J Cij 

= - dkiidki In det U {k) 
2 Jcij 

= Ntt, {N: integer ). (131) 

Thus the Berry phases e ' takes only two possible values, and we use it as the 
topological index, 

(_l)-[cd=e^fc.,^^^('=)^±l. (132) 

Like the one-dimensional winding number wid in §3.3, the topological index (^—lyl'^ij^ 
characterizes the "one-dimensional gapped superconductor" along the path Cij. When 
(^—lyi^ij] — —1, the bulk-edge correspondence yields that there exists a zero energy 
edge state at the corresponding momentum in the surface BZ. On the other hand, 
when (— 1)'^1'^»j1 = 1, no ZEABS is required at the corresponding momentum in the 
surface BZ. 

Up to now, we have not used the specific property of odd-parity superconductors. For 
odd-parity superconductors, the combination of the inversion and the electromagnetic 
U{1) gauge symmetry, Cka iPaa'Cka', is manifestly preserved. Thus, H{k) has the 
following symmetry 

U.^H(k)U = H(-k), n= ( ^ ° I. (133) 

\ -P* J 

From this, we have [i7(Fj), 11] = for the time- reversal invariant momentum Fj. Thus, 
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the quasiparticle state \un(ri)) at Fj is simultaneously an eigenstate of U, 

n|«„(r,)) = 7r„(r,)|«„(r,)). (i34) 

Then, the product of 7r„(rj) for occupied states, Y[n<o'^n(^i)^ following inter- 

esting properties: (a) it takes only discrete values rin<o'^"(^«) ~ (^) value 

can change only when the gap of the system closes at Fj. These properties are what we 
expect for the topological index, and suggest that nn<o^«(-'^*) ^ ^^T^d of topological 
index. Indeed, we can relate nn<o^^(-'^*) topological index (— Ij'^P^il as,^^^^ 

(_i)^[cd^ jj^„(r,)^^(r,). (135) 

n<0 

For ordinary superconductors, the superconducting gap is much smaller than the 
Fermi energy. Therefore, we reasonably assume that the typical energy scale of the gap 
function A(Fj) at the time-reversal invariant momentum is much smaller than that of 
£^(Fj). Under this assumption, we can deform the gap function A(fe) adiabatically so as 
A(Fj) — > without closing the bulk energy gap. Because of the topological nature of 
(^—lylCij]^ this adiabatic process does not change the value of (— 1)''[c<j1. in the process 
A(Fi) ^ 0, the BdG Hamiltonian at F^ reduces to //(F^) ^ diag(£(Fi), -£^(Fi)). For 
this simple Hamiltonian, we can evaluate 7r„(Fj) rather easily, and we obtain the final 
expression, ^^^^ 

(-l)-[Cd ^ sgn£,(F,)sgn£,(F,), (136) 

a 

where ea{^i) is an eigenvalue of i^(Fj) and the product of a is taken for all eigenstates of 
£(Fi). This formula implies that for odd-parity superconductors the gapless boundary 
state can be predicted by the Fermi surface structure. 

When the odd-parity superconductors are fully-gapped in two dimensions, the Chern 
number Ci also characterizes the topological phases. By using eqs. (48) and (128), it 
can be hnked to the topological index vlCij] as 

Ci = ^ /" d^kB,{k) = - [ d^kB,{k) 

^ -i dhA^r\k) = u[C^2]-i^[G^^l (137) 

where T| is the upper half of T^, and Cjj is the time-reversal invariant path illustrated 
in Fig. 34(b). Therefore, eq.(136) yields 

i-^f"- n Sgn£a(^^)■ (138) 

a,i=l,2,3,4 
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One can confirm the validity of tlie formulas (136) and (138) by applying them to 
the spinless superconductors in Fig. 27. For example, from the Fermi surface structures 
illustrated in the top row, we find that only the case of Fig. 27 (b) has an odd Chern 
number. Correspondingly, a single edge mode appears on each edge in Fig.27 (b). Fur- 
thermore, it is found that the difference between Fig.27 (a) and (c) can be understood 
as the difference in (-1)'^[C!.,] g^^^j^ ^^gg. Qne find that (-l)'^[Ci2] = (_i)''[C34] = _i 
for Fig. 27 (c), while ( — l)''Pi2] = (_ i)''[C34] _ ^ ^^j. (^^^^ ^^le zero energy edge states in 
(c) are originated from the non-trivial values of these topological indices. 



k2 




ise 



Fig. 34. (Color online) The time- reversal invariant momenta Fi, and the time- reversal invariant 
closed path passing through Ti and Tj in the BZ. a) ID BZ. The solid line denotes Ci2- b) 2D 
BZ T^. c) 3D BZ and the surface BZ of a (100) face, d) 2D band structure for a slab with a (100) 
face for the 3D time-reversal invariant odd parity superconductor that has the Fermi surface with 
(_1)£>[C34] = (_i)i>[C56] = _i and the gap function A(fe) = id{k) ■ aoy with d,(fc) = Asin/ci. The 
insets show the Fermi surface (top) and the 2D surface BZ (bottom). [Reproduced from Fig.l of Phys. 
Rev. B 81, 220504(R) (2010) by Sato] 

Now consider the time-reversal invariant case. The time-reversal invariance implies 
that (— l)'^[^'j] is always trivial, i.e. (— 1)'^[^»j] = 1. However, use of the time-reversal 
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invariance as well as the particle-hole symmetry makes it possible to define another 
topological index. 

Because of the time-reversal invariance, the quasiparticle states form Kramers pairs, 
I'^ni^)) ~ I) II)) \'^n{f^)) = T\ull{—k)) , with the time-reversal operator T. The new 
topological index is defined by the gauge field for the "half" of the Kramers doublets, 
say, Al^~\k) = i X^„<o('"n(^)l^feil'"n('^))- Then the topological index (^—lyi'^ij] is defined 
as 

(_l)^[Cd^e'^«^.-^^'"'^'=\ (139) 

In a manner similar to the above, for odd-parity superconductors, (— are deter- 
mined by the Fermi surface structure,^^^^ 

(-lyidA ^ sgn£2a(r,)sgn£2a(r,), (140) 

a 

where £a(ri) is an eigenvalue of £{ri), and we have set £20(1^1) = £2a+i(ri) by using 
the Kramers degeneracy. 

Furthermore, other topological indices characterizing time-reversal invariant su- 
perconductors can be finked to the Fermi surface topology. For the full-gapped two- 
dimensional case, their topological nature is characterized by the Z2 number (^—ly^dTi 
introduced in §3.5 and §3.6. For odd-parity superconductors, it satisfies^^^'^^^^ 

(-l)-2dTi^ n Sgn£2a(rO, (141) 
a,i=l,2,3,4 

where is the time-reversal invariant momenta in Fig.34(b). Thus, we can predict the 
helical Majorana edge mode from the knowledge of the Fermi surface. In addition, for 
the full-gapped three-dimensional case, we have another topological index known as the 
three-dimensional winding number u^sa.^^^'^^^'^^^'^^^^ (When w^^ takes a nonzero integer 
n, there exist n Majorana cones on the surface.) We can also connect this to the Fermi 
surface structure as^^^'^^^'^^^) 

(-l)-3d= Yl sgn£2a(ri), (142) 

a,i=l,--- ,8 

with Fj in Fig. 34 (c). The last formula determines the parity of the three-dimensional 
winding number, thus it gives a sufficient condition for nonzero Wsa- In Fig.34 (d), 
we illustrate surface states of a three-dimensional time-reversal invariant topological 
superconductor. The Fermi surface of this model leads to (— l)"'^'! — 1, but (— l)^[C34]=-i 
and (— l)'^[C56l=-i Consistently, we have two Majorana cones on the surface at the time- 
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reversal invariant points and F^q in the surface Brillouin zone. (For the definition 
of Fij, see Fig 34(c).) 

We would like to emphasis here that the correspondence between the Fermi surface 
and the gaplcss surface state discussed in the above are inherent to spin-triplet/odd- 
parity superconductors. Thus we may identify the spin-triplet/odd-parity superconduc- 
tivity through the direct measurement of the surface state, irrespective of the details of 
the gap function. For instance, the formula (142) has been applied to Cuj;Bi2Se3 
in order to identify its possible odd-parity superconductivity.^^^^ Detailed information 
on ABS for such a material can be predicted by using the formulas (140) and (141) as 
well. 

3.8 Majorana fermion on the surface of topological insulator 

As discussed in §3.5, the topological insulator is a new state of matter with time- 
rcvcrsal (T— symmetry). ^^^"^^^'^^^^ Especially, on the surface of the 3D strong topological 
insulator (STl), there appears odd number of 2D Dirac fcrmions. Even with a disorder 
scattering etc., at least one Dirac fermion is protected by the gap and topology of 
the bulk states and is stable as long as the perturbation is T— symmetric and does 
not destroy the bulk gap. In the case of Bi2Se3 and Bi2Te3,^^^^ there appears only one 
surface Dirac fermion, which can be regarded as the "half of the 2D electrons because 
the spin direction is determined by the momentum. This " fractionalization" of the 
electrons is one of the key concepts in the topologically ordered states. This situation 
is similar to the case of spinless fermion, and when the Cooper paring occurs, there is a 
possibility that Majorana fermions will appear. ^^^^ Fu and Kane were the first to study 
the superconductivity of the surface Dirac fcrmions of the STl induced by the proximity 
effect by the superconductors attached to it.^^*^) Also they considered the ferromagnetic 
insulator (FI) put on the STl, and its interface with the conventional superconductor 
(S). They predicted the emergence of the one-dimensional chiral Majorana mode as 
an Andreev bound state.^^^^ This situation is similar to that discussed in §3.4. Also 
by applying the magnetic field, the vortices penetrate the sample. At the core of each 
vortex, the Majorana bound state is reahzed at zero energy. ^^^'^^^^ 

Majorana interferometry has been proposed to detect the Majorana fermions. 
The schematic configuration is shown in Fig. 35 . Everything is on top of the STl, and 
there appears the chiral fermion channel at the interface of two magnetic domains with 
t and 4- magnetization, while the chiral Majorana fermions appear at the interface 
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Fig. 35. Majorana interferometer on top of a strong topological insulator. [Reproduced from Fig. 1 
of Phys. Rev. Lett. 102 216403 (2009) by Fu and Kane. ] 



between FI and S. Let and c being the creation and annihilation operators of this 
chiral fermion. At the branching points, this chiral fermion is connected to the two 
Majorana edge channels for upper and lower semicircles. Let 71 (72) be the operator of 
the Majorana fermion in the upper (lower) semicircle. Let us consider the process that 
an electron is entering from the left to the branching point described by 

c-^ 71 + 272. (143) 

At the right branching point, these two Majorana fermions merge into an electron as 

7b 72 c = 71 + 272 (144) 

or 

7b 72 c"^ = 7i - ^72 (145) 

depending on the even or odd number of vortices $/0o (0o = hc/2e) trapped by the 
superconducting donut according to the Aharanov-Bohm effect. This leads to strong 
dependence of the conductance G between the left and right leads. Namely, G = 2e^ /h 
when $ = 2n X 00 while G = when $ = {2n + 1) x 0o with n being an integer. This 
prediction is based on the fact that the Majorana fermion contains the particle and hole 
components with equal weight, or an electron is composed of two Majorana fermions. 

The other direction of the research is to study the possible manipulation of the Ma- 
jorana fermion and its relevance to the transport properties in the transverse direction 
to the interface, For that purpose we consider the Hamiltonian for the surface 
state on STI influenced by the superconductor and ferromagnet ( Fig. 36 ) as given by 

( H{k) + M A \ ^ ^ 

Hs=\ ^ \ . ^ . ] 146 

I -A* -H* {-k) - M* j 

where H{k) = vpio'xkx + o'yky) — ^[<d{~x) + <d{x — d)] describes the 2D Dirac fermion. 
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and M = m ■ &Q{d — x)Q{x) is the magnetic exchange interaction with m & = mxd'x + 
mydy + mzCTz- Here, /i, a, vp, m denote chemical potential, Pauli matrices, velocity, 
and magnetization (times the exchange coupling constant which we assume to be 1), 
respectively. Interestingly, m couples to the Dirac fermion as an effective vector 
potential A of the electromagnetic field. For the moment, we assume that the pairing 
symmetry of superconductor is spin-singlet s-wave and A is given as A = i&yAQlx — d) 
and A = iay[Ae{x - d) + Ae(-x) exp(zy?)] for N/FI/S (Fig. 36 (a)) and S/FI/S 
(Fig. 36 (b)) junctions formed on the surface of STI, respectively, where ip denotes 
the macroscopic phase difference between left and right superconductor. (N represents 
the normal metal which is the 2D Dirac fermion in the present case.) We also assume 
that the magnitude of A is smaller than that of the bulk energy gap of superconductor 
deposited on STI because the S/TI interface is not ideal.^^^^ 



(a) z (b) - m 




Fig. 36. (Color online) Schematic illustration of the junction, (a) Normal metal (N) / Ferromagnetic 
insulator (FI) / Superconductor (S) junction and (b) S/FI/S junction formed on the surface of 3D 
topological insulator (STI). The current is flowing on the surface of STI. [Reproduced from Fig. 1 of 
Phys. Rev. Lett. 103 107002 (2009) by Tanaka et. al] 

First, we solved the Andreev bound state (ABS) at the interface of FI/S. As an exam- 
ple, we show in Fig. 37 the energy dispersion as a function of the angle 6 (corresponding 
to the momentum along the interface) in the right panel, and the corresponding con- 
ductance a of the N/FI/S junction in the left panel for several values of fi/rriz. Note 
that the slope of the energy dispersion affects the conductance a sensitively. Therefore, 
it is concluded that the ID Majorana channel can be controlled by the magnetization 
of the FI and the chemical potential /i. 

Next, we consider the Josephson effect in S/FI/S junction (Fig. 36(b)). This is 
mainly due to the tunneling process between the two Majorana edge channels at the 
two interfaces between S and FI. One interesting prediction is that the current-phase 
relation can be shifted by the rrix. Figure 38 shows the Josephson current as a function 
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Fig. 37. (Color online) Left Panel: Chiral Majorana mode energy dispersion Ei, as a function of the 
incident angle 0. Right Panel: Normalized tunneling conductance a in N/FI/S junctions, vrizd/vp = 1 
and my/mz = 0. a: ^/friz — 1, b: ^/rriz — 2 and c: fi/rriz — 0.5. [Reproduced from Fig. 2 of Phys. Rev. 
Lett. 103 107002 (2009) by Tanaka et. al.] 

of the phase different (p between the two superconductors for different values of m^, 
where the parallel shift of the curves is seen. 

Josephson Current 




Fig. 38. (Color online) Josephson current in S/FI/S junctions is plotted with a: rrix/mz = 0, b: 
f^xI'Tfiz = 0.4 and c: m^lrriz = —0.4. We choose nizd/vp = 1, [ijmz = 1 and my/mz = 0. T = 0.05Tc 
with transition temperature Tc- [Reproduced from Fig. 4 from Phys. Rev. Lett. 103 107002 (2009) by 
Tanaka et. al.] 

The generalization to the anisotropic pairing superconductors on STI has been dis- 
cussed in i-ef.225,226) _ jg found that the proximity effect to the superconductivity is 
strongly dependent on the pairing symmetry. For spin-triplet p-wave superconductor, 
there occurs no proximity effect because the spin direction of the 2D Dirac electrons 
is not consistent with that of the p-wave pairing. In the cases of spin-singlet s-wave 
and d-wave pairings, there occurs the proximity effect, and the Majorana fermions ap- 
pear at the interface of FI and S. For spin-singlet dxy-pa.inng, the zero energy states 
appear along the edge of the sample even without the magnetization of FI, i.e., with- 
out T— symmetry breaking, which has been well-known. However, the new aspect here 
is that the zero energy dispersionless bound states are Majorana fermions due to the 
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Table IX. Edge states of various superconductors. Each column describes the type of bulk supercon- 
ductor, zero energy Andreev bound state (ZEABS), the presence or absence of the energy dispersion, 
topological index characterizing the ZEABS, symmetry of the induced odd-frequency pairing, the type 
of Majorana fermion, and the corresponding sections in the text, respectively. 

proximity to STL With the 2;-component of the magnetization, the Majorana fermion 
obtains the dispersion, which results in the very different behavior of the conductance 
compared with the case of spin-singlet s-wave superconductors. 

4. Discussion 

4-1 Interplay between symmetry and topology 

Up to now, we have studied induced odd-frequency pairings in §2 and topological 
bulk-edge correspondence in §3. In this section, we discuss the relation between these 
two. 

First, all the known zero energy edge states at the surface of superconductors can 
be characterized by the topological indices, and are topologically protected. These zero 
energy states are always associated with the pure odd-frequency pairing at the surface. 
(When the ABS has the energy dispersion, this statement apphes only to the sector of 
ky — 0.) In Table IX, we summarize the results for various situations studied thus far. 

For s-wave superconductor, which belongs to ESE, there is no ZEABS because 
it is topologically trivial. At the same time, odd-frequency pairing is absent at the 
surface. However, dxy-wave superconductor, which also belongs to ESE, has ZEABS 
characterized by the topological index wid- wid is the winding number defined for the 
one-dimensional systems when ky is fixed as discussed in §3.3. The formation of the 
ZEABS stems from the sign change of the pair potential felt by injected and reflected 
quasiparticles. The similar situation occurs for p-^-wave superconductor which belongs 
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to ETO pairing. In these cases, pure odd-frequency pair amplitudes (OSO for dxy, OTE 
for exist at the surface. The ZEABSs can be regarded as double Majorana modes 
satisfying 

7.(y=7-a(-y, (<T=t4)- (147) 

Next, we discuss the spinless chiral p-wave superconductor where the time-reversal 
symmetry is broken in the bulk. There exists ABS, however, in contrast to the cases 
above, it has a dispersion linear in ky. Because of the broken time-reversal symmetry, it is 
characterized by the Chcrn number W2d in a manner similar to the quantum Hall states. 
The induced odd-frequency pair amplitude at the surface belongs to OTE symmetry. 
The present ABS is a chiral Majorana mode, which is a genuine Majorana fermion. The 
case of the spinful chiral p-wave superconductor can be regarded as the two copies of the 
spinless case discussed above. Lastly, the noncentrosymmetric superconductors (NCSs) 
are discussed. In these superconductors, the parity of the pair amplitudes are mixed. 
For the case of s-wave and p-wave mixture, there are two topologically distinct states 
depending on the relative magnitudes of the corresponding pair potentials. When p-wave 
component is dominant, the ABS is generated as a helical Majorana edge mode. The 
pair amplitude induced at the surface has OTE+OSO symmetry. The topological index 
characterizing the ABS is the Z2 invariant. On the other hand, when s-wave component 
is dominant, no ABS is generated because it is smoothly connected to the topologically 
trivial s-wave superconductor case without closing the bulk gap. For {dxy + p)-wave 
pair potential case, a single branch of dispersionless ABS appears if ky satisfies some 
condition. This ZEABS is a single Majorana fermion. The topological index W defined 
by eq.(123) ensures its topological stability. ^^^^ The symmetry of induced odd- frequency 
pair amplitude is OSO-^OTE. 

4-2 Related topics 

There are several important issues which are not addressed in this review article 
due to the limited space. We will briefly mention each of them by indicating relevant 
references. 

In this review article, we mainly focus on the interface and surface of supercon- 
ductors. There are several works on the bulk superconductors from the viewpoint of 
topology and/or odd-frequency pairings. Non-interacting topological insulators and 
topological superconductors have been classified based on the universal class of the 
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Anderson localization. ^^^'^^'''^ The similar results have been reached independently by 
Kitaev^^^) in the framework of K-theory. The classification scheme also has been apphed 
for solitons in topological insulators and topological superconductors. ^^^^ Possible bulk 
odd-frequency superconductors have been an important issue in many-body physics. 
Recently, there are several proposals of odd-frequency gap functions in quasi one- 
dimensional systems/^' magnetic superconductors, heavy fcrmion systems, 
and strongly coupled electron-phonon systems. However, it is noted that the thermo- 
dynamic stability and physical properties of the bulk odd-frequency pairing supercon- 
ductors are still controversial issues.^^'^^-*^'^^^) 

Since superfluid ^He is a well established spin-triplet p-wave superfluid, it is natural 
to expect ABS.^^'^^) ABSs in superfluid ^He have been predicted in the 1980's.^^'^^) 
However, the detection of them is not easy since Cooper pairs of superfluid ^He do not 
have electric charge. Surface ABS of ^He B phase has been observed experimentally 
by the transverse acoustic impcdancc^^'^^ in 2005 and specific heat mcasurcmcnts^^^^ in 
2006. In the actual surface of ^He B phase, diffusive scattering by surface roughness 
plays an important role.^^^'^^^^ When a surface is specular, SABS has a Dirac like hnear 
dispersion which have been recently interpreted as Majorana fermion.^^^'^^^^ There are 
several theoretical discussions about SABS from the viewpoint of topological supercon- 
ductivity and Majorana fermions.^^^"^^^'^^^"^^^^ It is noted that surface Majorana cone 
has been experimentally observed for highly specular wall of ^He B phase. ^^^'^^^^ Sur- 
face Majorana ABS (SMABS) has been also discussed in slab geometry of ^He A and 
B phases. ^^^'^"^^^ An analogous system of diffusive normal metal/spin-triplet supercon- 
ductor can be realized in superfiuid aerogel/'^Hc hybrid system. Anomalous proximity 
effects discussed in §2.4 are also expected in the aerogel, since it can be regarded as a 
diffusive Fermi liquid. ^^^^ 

Recently, promising candidate of topological superconductor has been discovered 
in doped topological insulator Cuj;Bi2Se3.^^^^ Several experimentaP^^' ^^^^ and theoreti- 
j,g^pi4,249) researches have started. MSABS are also expected in the present new super- 
conductors. 

Majorana fermion can be a platform of topological quantum computation. Var- 
ious systems have been proposed to produce Majorana fermion; spin-triplet p-wave 
superconductors / supcrfiuids,^'''^'^'''^'^'''^'^^^'^^^"^^^^ z/ = 5/2 fractional quantum Hall sys- 
tem, ^^^^ spin singlet s-wave superconductor on topological insulator surface with vor- 
tex or the interface with ferromagnet,^^^'^^^'^^^^ noncentrosymmetric superconductors 
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with Zeeman field/^^'^^^'^^^^ s-wave superfluid of cold atoms with laser generated spin- 
orbit interaction/^^' semiconductor/SC systems including Rashba spin-orbit cou- 
pling with FM-'^^''') or Zeeman field, one-dimensional system made of 2D QSHS^^°) 

or nanowire,^^^"^^^^ quantum anomalous Hall/ SC hybrid system,^^^'^^^) and Half metal 
/SC system.268.269) 

The detection of Majorana fcrmion is not a trivial issue since it is a neutral par- 
ticle. Majorana fermions in condensed matter physics are emergent particles result- 
ing from the many-body interactions of electrons which are charged. Thus, they can 
be detected through the conversion processes to electron influenced by the phases of 
the superconducting order parameters. In addition to the Majorana interferometry on 
topological insulator discussed in §3.8,^^°"^^^^ there are many other proposals for this 
purpose: (i) Josephson effect,^^°'^^^'^^^'^^°"^^^^ (ii) charge transport thorough Majorana 
fermion,^^^'^'''^"^''"'') (iii) transport noise,^''®'^™^ and (iv) spin anisotropy.^^^'^°^'^^''''^^^'^^°) 

4-3 Perspectives 

The two major fields opened up in 80's are the quantum Hall effect and the high- 
Tc superconductors. The former introduced the concept of topology and the latter the 
electron correlation in condensed matter physics. The idea of the topology extended 
to the time-reversal symmetric systems and the topological insulators have been dis- 
covered. On the other hand, the electron correlation naturally leads to the anisotropic 
Cooper pairings with various symmetries. These two streams are now merging into one 
producing many interesting new directions in the current condensed matter physics. 

It is now recognized that there are 3 distinct classes of currents in solids; (i) One is 
the conventional Ohmic currents with dissipation, (ii) The second is the superconduct- 
ing currents and/or the superfluid flow with broken gauge symmetry and macroscopic 
quantum coherence, (iii) The third one is the "topological currents" induced by the 
nontrivial topology of wavef unctions, and does not require the broken gauge symmetry. 
The last two have the common feature of dissipationless currents. The last one has not 
been paid attention thus far, but will be very unique because it can operate at room 
temperature when the gap protecting the topology is larger than its energy scale. 

These three classes have been considered to be rather clearly separated. However, 
as we have described in this review article, class (ii) and class (iii) are overlapping each 
other. Namely, there is a class of topological supercurrent at the edge/surface of the 
superconductors/superfluid. Therefore, the interplay of the broken gauge symmetry and 
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topology offers an intriguing new arena for the future condensed matter physics. 

These fundamental understandings of the currents in solids will be essential to de- 
velop the electronics in the next generation. The dissipationless operations of devices 
are the key ingredient for the future applications to electronics, and hopefully the topo- 
logical supercurrent will serve this goal. The other important issue is the quantum 
computing. Majorana fermion is regarded as the promising candidate for fault-tolerant 
qubit. In general, the topological orders bring about the electron fractionalization, which 
gives the robustness of the quantum states. The three degrees of freedom of electrons, 
i.e., spin, chirality, particle-antiparticle, can be fractionalized at the edge/surface of 
topologically ordered bulk states, and chiral Majorana fermions can be regarded as the 
ultimate fractionalized quantum state. On the other hand, the other direction is to 
recombine these fractionalized particles into new types of particles. ^^^^ The more the 
degrees of freedom, the more the system is susceptible to the external stimuli, and hence 
we can adjust the robustness and sensitivity of the system at will by manipulating the 
fractionalization / recombination. 

To verify these theoretical proposals, experimental progresses are indispensable. 
There are several relevant experimental techniques such as //SR, NMR, neutron scat- 
tering. X-ray diffraction, and magneto transport. In particular, the direct observation 
of edge/surface states by the tunnehng spectroscopy and ARPES is highly desired. For 
more details, the readers are referred to other contributions in this special volume. 

Finally, we would like to stress that physics of interface/surface superconductiv- 
ity includes unexplored quantum phenomena and will be a central topic in condensed 
matter physics in the decades to come. 
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